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ABSTRACT 


This  work  is  concerned  with  predicting  the  forces  acting  on  slender 
bodies,  namely  hydrofoils  and  wedges,  in  rotational,  super cavi tat ing 
flow.  Methods  are  given  for  establishing  not  only  qvialitative  but 
quantitative  measures  of  the  effects  of  rotation  in  linearized,  super- 
cavitating  flows. 

A  linearized  theory  is  developed  for  steady,  two-dimensional 
flow  under  the  assumption  that  the  flow  has  a  constant  vorticity 
throughout.  The  effects  of  gravity,  viscosity,  and  surface  tension 
are  neglected.  Tialin's  original  closed-cavity  model  is  employed.  A 
basic  assumption  of  the  theory  is  that  the  slender  body-cavity  combina¬ 
tion  caxises  only  small  perttarbations  in  the  velocity  components  of 
the  basic  shear  flow.  The  stream  function  of  the  rotational  flow 
satisfies  Poisson's  eqaiatlon,  which  is  a  linear,  inhomogeneovis, 
partial  differential  eqaaation.  By  using  a  particular  solution  of  this 
equation,  the  linearized,  rotational  problem  is  reduced  to  a  problem 
Involving  Laplace's  equation  and  harmonic  perturbation  velocities. 

The  bomdary  conditions  for  the  perturbation  velocities  are  established 
from  facts  known  about  the  body-cavity  combination  in  the  supercavitat- 
ing  shear  flow.  The  resulting  boundary  value  problem  is  solved  by 
the  use  of  conformal  mapping  and  singularities  from  thin  airfoil 
theory. 

The  theory  is  applied  to  asymmetric  shear  flow  past  wedges  and 
hydrofoils  and  to  symmetric  shear  flow  past  wedges.  Analytic  expressions 
are  given  for  pressure,  drag,  lift,  and  moment  coefficients  as  well  as 
cavity  length,  cavity  area,  and  cavitation  nuniber  relationships.  The 
presence  of  vorticity  is  shown  to  create  significant  changes  in  those 
forces  acting  on  the  slender  bodies  and  in  the  shape  and  size  of  the 
trailing  cavities.  The  results  are  summarized  in  tables,  graphs,  and 
tabulated  numerical  calculations. 
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1.  IMTRODUCTION 


Cavitation  occurs  in  a  fluid  flow  as  a  consequence  of  local  pressure 
reduction,  generally  brou^t  about  by  high  local  velocities.  The  devel¬ 
opment  of  high-speed  submarines,  underwater  missiles,  and  other  vehicles, 
together  with  the  surface-piercing  hydrofoil  ship,  has  renewed  interest 
in  the  large  scale  effects  of  cavitation.  The  hydrofoil  and  wedge  (or 
strut)  are  practical  parts  of  the  total  hydrodynamic  system  of  most  of 
these  vehicles.  In  many  cases  these  parts  have  long  and  slender  cross 
sections  with  their  greatest  dimension  nearly  parallel  to  the  flow  direc¬ 
tion.  At  sufficiently  high  speeds,  common  for  present  vehicles,  these 
slender  bodies  produce  long,  trailing,  steady-state  cavities  as  the  result 
of  air  ventilation  or  cavitation.  The  characteristics  of  these  so-called 
supercavitating  flows  about  bodies  are  of  particular  interest  to  the  de¬ 
sign  engineer.  If  the  cavity  pressure  does  not  differ  greatly  from  the 
free-stream  static  pressure,  the  velocities  near  the  body  and  cavity  do 
not  differ  greatly  from  the  free-stream  speed.  It  is  possible,  then,  to 
study  the  flow  by  means  of  a  linearized  theory  which  is  based  on  the  well- 
known  two-dimensional  theory  of  thin  airfoils. 

Tulin  [1]^  appears  to  be  the  first  to  have  used  the  linearized  theory. 
The  work  published  since  the  appearance  of  Tulin's  paper  has  been  concerned 
with  both  higher  order  linearized  theory  [2]  and  extensions  of  the  first 
order  theory  to  include  effects  found  in  real  flows,  e.g.,  surface  tension 
[3]  and  gravity  [4,5].  The  linearized  theory  has  been  applied  to  many 
practical  problems  which  were  insoluble  by  more  classical  means.  The 
progress  of  this  work  up  to  i960  is  summarized  in  three  papers,  two  by 
Tulin  [6,7]  and  one  by  Parkin  [8]. 

In  the  linearized,  two-dimensional  theory  the  effects  of  viscosity 
are  usually  neglected.  The  flow  is  assumed  to  be  irrotational  and  the 
velocity  is  assumed  to  be  uniform  for  points  far  from  the  slender  body. 
However,  since  no  fluid  is  completely  Inviocid  all  real  flows  are  rota¬ 
tional.  Even  when-  viscosity  is  neglected,  the  flow  picture  may  sometimes 
be  best  represented  by  a  rotationsil  flow. 

4!he  numbers  in  brackets  refer  to  the  references  listed  at  the  end  of  the 
work. 
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Many  rotational  flows  have  already  been  studied  empirically  and  ana¬ 
lytically.  For  example,  the  only  known  exact  solution  for  the  problem  of 
finite  wave  motion  is  Gerstner's  trocholdal  wave  which,  while  producing 
a  rotational  flow,  also  satisfies  exactly  the  constant  pressure  boxindary 
condition  at  the  free  surface.  Also,  two  common  physical  flows,  the 
eye  of  a  typhoon  (a  forced  vortex)  and  uniform  viscous  open  channel  flow, 
are  rotational.  Lamb  [9]  Groen  [10]  have  studied  another  rotational 
flow  -  the  propagation  of  small  surface  waves  on  a  stratified  fluid.  The 
rotational  flows  studied  by  Yih  [11]  are  of  partlculso'  Interest  because 
his  results  show  the  Importance  of  vorticity  in  reproducing  physical  ef¬ 
fects.  He  considers  the  steady,  rotational  flow  of  an  inviscld  fluid  in 
a  two-dimensional  channel  or  a  circular  tube  toweurd  a  sink.  His  solutions 
show  the  imusiial  (for  inviscld  theory)  features  of  separating  streamlines 
and  comer  eddies.  Note  that  for  a  viscous  fluid  flowing  in  a  channel  or 
pipe  with  an  abrupt  contraction,  eddies  occur  in  the  comers  formed  by 
that  contraction,  but  such  eddies  are  not  predicted  by  an  irrotational 
analysis.  In  19^3>  Tsien  [12]  recognized  that  there  were  many  applications 
in  two-dimensional  airfoil  theory  \diere  irrotational  flow  conditions  are 
not  satisfied.  He  states  in  his  paper  on  airfoils  in  shear  flow  that,  for 
example,  the  large  vertical  velocity  gradient  near  the  ground  can  be  ap¬ 
proximated  to  the  first  order  by  a  flow  with  a  linear  velocity  distribu¬ 
tion  (a  shear  flow).  Thus,  according  to  Kronauer,  "...  the  discussion 
suggests  that  over  a  limited  stream  length  the  essential  character  of  the 
motion  may  be  closely  approximated  by  specifically  neglecting  the  viscous 
forces  acting  in  that  stream  length,  but  by  Inclxidlng  (pexhaps  approxi¬ 
mately)  the  effects  of  viscotis  forces  up-stream  [13  ] . " 

The  present  work  may  be  regarded  as  an  extension  of  both  Tulin's 
llneeurlzed  theory  for  supercavltatlng  flow  and  Tsien 's  method  for  rotational, 
non-cavltatlng  flow.  Methods  are  given  for  establishing  not  only  qualita¬ 
tive  but  quantitative  measures  of  the  effects  of  rotation  in  linearized, 
supercavltatlng  flows.  The  effects  of  gravity  and  surface  tension  are  spe¬ 
cifically  neglected.  In  the  study  a  uniform,  parallel,  Irrotational  flow 
is  perturbed  by  a  uniform  shear  flow  -  the  simplest  perturbation  of  the 
parallel  flow  and  a  flow  with  constant  vorticity.  The  vorticity  is  pre¬ 
sumed  to  have  come  from  some  up-stream  disturbance,  e.g.,  a  boundary  layer 
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developing  on  the  body  of  a  vessel  or  test  tunnel  walls.  Hydrofoils  and 
struts  which  lie  in  the  slipstream  or  wake  of  other  con5»onent8  woxild  be 
in  such  a  vortex  field. 

The  perturbed  flow  and  the  physical  problems  to  be  studied  are  shown  in 
Figure  1.  The  perturbed  flow  is  also  a  \iniform  shear  flow  and,  as  such, 
is  characterized  by  a  linear  velocity  distribution  and  a  constant  vortl- 
city  e  throughout.  The  irrotational  and  rotational  perturbed  flows  both 
satisfy  the  equation  of  continuity;  therefore  the  stream  function  ^  exists 
in  both.  However,  the  velocity  potential  0  can,  of  course,  exist  only 
in  the  irrotational  flow.  Physically,  the  flow  studied  is  two-dimensional, 
and  the  fluid  is  of  infinite  extent.  The  flow  is  supposed  to  detach  at 
the  leading  and  trailing  edges  of  the  hydrofoil  and  at  both  edges  of  the 
blunt  base  of  the  symmetric  wedge.  Finally,  it  is  assumed  that  the  cavity 
length,  when  measured  from  the  leading  edge  of  the  body,  is  greater  than 
that  of  the  solid  body,  i.e.,  full  cavitation  occurs. 

The  main  paper  considers  rotational,  supercavitating  flow  past  wedges 
and  hydrofoils.  The  bodies  chosen  have  been  restricted  to  the  slender 
wedge  and  the  flat-plate  hydrofoil.  The  Appendixes  include  a  discussion 
of  supercavitating  flow  past  wedges  in  a  transverse  gravity  field.  The 
solution  to  this  problem  complements  the  gravity  flow  solutions  already 
given  by  Parking  [4]  and  Acosta  [5]  and  arises  directly  from  the  solution 
of  the  rotational  flow  past  a  wedge. 
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2.  THE  LIHEARIZEn}  THEORY 


In  this  section,  a  linearized  theory  is  developed  for  two-dimensional, 
supercavitating  flows  under  the  assumption  that  the  flow  has  a  constant 
vorticity  throvighout.  The  effects  of  gravity  and  surface  tension  are 
neglected.  The  flows  considered  are  those  past  slender  wedges  and  flat- 
plate  hydrofoils  as  shown  in  rigure  1.  It  is  further  assumed  that  the  flow 
is  steady  and  that  the  fluid  is  both  incompressible  and  inviscid.  Because 
of  its  simplicity  and  convenience,  Tulin's  original  closed-cavity  model 
is  employed  in  this  work.  It  should  be  noted  that  Other  linearized  models 
are  available  [l4, 15]. 

2.1,  Notation  and  Boundary  Conditions 

Figure  2  shows  two  typical  fully  cavitated  flows;  the  notation  vised 
is  introduced  in  the  following  discussion.  The  incompressible  fluid  has 
a  constant  density  p  .  The  base  flow  which  is  a  parallel,  uniform  shear 
flow  has  been  distvirbed  by  the  introduction  of  the  slender  body  of  unit 
length.  The  vinit  length  body  is  used  without  loss  of  generality,  since 
it  is  equivalent  to  normalizing  the  problem  on  the  actual  body  length. 
Although  it  is  certainly  disturbed  in  the  nei^borhood  of  the  slender 
body,  the  base  flow  is  assumed  to  be  undisturbed  at  infinity  and  is  char¬ 
acterized  there  by  a  constant  vorticity  e.  The  wedges  are  aligned  symmet¬ 
rically  with  their  longitudinal  center-lines  parallel  to  the  x-axis,  and 
hydrofoils  are  placed  at  an  angle  of  attack  a  with  respect  to  the  x-axis. 

The  origin  of  the  rectangular  coordinates  is  at  the  leading  edge  of 
the  solid  body.  In  terms  of  these  coordinates,  the  velocity  profile  at 
X  =  -“  is  -  €  y,  with  U„  representing  the  velocity  at  ( -«,  o) .  The 
pressure  at  infinity  is  taken  to  be  the  undisturbed  static  pressure  1:^0. 

The  flow  velocities  U  and  V  are  in  the  x-  and  y-directions  respectively. 
The  total  velocity  at  any  point  in  the  fluid  is  q,  while  the  velocity  on 
the  cavity  surfaces  is  q^.  The  closed,  trailing  cavity  \diich  springs  from 
the  solid  body  is  characterized  by  a  length  t  -  greater  than  one,  an  ordi¬ 
nate  y  (x),  and  a  uniform,  constant  pressure  p  .  It  is  presumed  that 
c  c 

the  cavity  is  filled  with  air  or  water  vapor.  The  pressure  p^  is  always 
less  than  or  equal  to  p„. 


I 
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In  general,  the  only  restriction  imposed  on  the  shape  yQ(x)  of  the 
slender  bodies  is  that  the  flow  over  the  body  must  satisfy  a  Brillouin- 
Villat  separation  condition  [l6 ] .  Under  this  condition  the  maximum  velo¬ 
city  on  the  surface  of  the  body  must  occur  at  the  separation  point;  this 
corresponds  to  requiring  a  fixed  and  Icnown  point  of  separation.  In  the 
flows  considered,  this  condition  is  satisfied. 

In  a  two-dimensional  flow,  the  vorticity  component  is 

5x  ^ 

and  the  rotation  oi  =  ^/2.  It  is  known  that  the  flow  has  a  constant  vor¬ 
ticity  ^  =  €  at  (-00,  y).  From  the  Helmholtz  theorem  on  the  permanence 
of  rotation,  it  follows  that  the  vorticity  €  persists  throvj^out  the 
fluid;  thus. 


^  by 


(2.1) 


Since  the  fluid  is  incompressible  and  the  flow  is  steady,  the  continuity 
equation 


bx  by 


=  0 


(2.2) 


must  be  satisfied  throughout  the  fluid.  In  a  rotational  flow  the  stream 
fxinction  =  \lf  (x,y)  exists  and  may  be  defined  so  that 


U  = 


b^ 


(2.3) 


The  function  \|f  satisfies  equation  (2.2)  identically,  and  Poisson's 
equation 


^  ^  .  7^*(x,y)  .  € 

Jbc  by 


(2.U) 
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l8  produced  When  \|f  Is  Introduced  into  equation  (2.1). 

The  Bernoulli  equation  for  rotational  flow  Is  foxmd  hy 
Integration  of  the  Euler  equations.  They  are.  In  this  case, 


U 


au 


U 


Sv 

3x 


oy 


Following  rearrangement,  these  eqxiatlons  "become 


Integration  gives  the  Bernoulli  equation  for  stesuiy  two-dimens ioneO. 
flow  with  constant  vorticity  (or  rotation), 


P  + 


pe^f  =  constant. 


(2.5) 


Poisson's  equation  (2,1+)  is  a  linear,  inhomogeneous,  partial  dif¬ 
ferential  equation.  Therefore,  by  the  use  of  superposition,  its 
solution  can  be  written  as  \jf  =  \lrg  +  )|fp,  with  =  0  and  V^i|fp  =  €  . 

The  stream  fi+nction  \lr„  represents  a  new  harmonic  flow;  i|f  is  a 

n  r 

particular  solution  of  equation  (2.1+),  From  equation  (2.3),  then, 

U(x,y)  =  Ug(x,y)  +  Up(x,y) 

V(x,y)  =  Vg(x,y)  +  Vp(x,y) 

If  one  lets  ii^p  =  =  €  as  desired,  and 

Up  =  -  6y, 


(2.6) 


I 


6 


u  «  Ug  -  €y 


Thus, 


It  Is  evident  that  any  problem  which  requires  a  solution  for  U  and  V 
can  be  reduced  to  an  equivalent  problem  for  the  harmonic  velocity 
and  V_. 

II 

With  the  flow  conditions  at  infinity  known,  it  remains  to 
establish  the  bovuadary  conditions  on  the  solid  body  and  cavity  walls. 
These  conditions  and  a  compatible  boundary  value  problem  are  estab¬ 
lished  in  terms  of  a  linearized  theory.  The  basic  assumption  of  such 
a  theory  is  that  a  slender  body-cavity  system  causes  only  small  per¬ 
turbations  in  the  velocity  coaqponents  of  the  basic  shear  flow.  This 
assumption  does  break  down  in  the  nele^borhood  of  singular  points. 

But,  as  Tulin  [6]  notes,  the  effect  of  this  breakdown  is  usually 
restricted  to  the  area  of  the  singular  points,  and  the  overall  char¬ 
acteristics  of  the  flow  are  often  well  reproduced  by  linearized 
solutions.  The  flow  pattera  is  a  conibinatlon  of  parallel,  uniform 
shear  flow  and,  superposed  on  this,  small  velocity  perturbations. 

Since  it  is  desired  to  formulate  eventually  an  harmonic  boundary  value 
problem,  the  flow  velocities  are  written  in  terms  of  the  harmonic  per¬ 
turbation  velocities  (u,v)  in  the  x-  and  y-  directions  respectively. 
From  equation  (2.6) 


components 


U  =  U  +  u  -  ey 

m 

V  e  V 

in  the  linearized  flow. 

In  a  linearized  theory  for  irrotatlonal  supercavitating  flow,  it 
is  assumed  that  the  perturbation  velocities  (u,  v),  the  attack  (or 
wedge  semi-)  angle  a,  the  body  shape  (or  camber) 

cavity  ordinate  are  small.  First  order  terms  in  these  quanti¬ 

ties  are  retained,  but  second  and  high  order  terms  are  neglected.  The 


► 


(2.6a) 
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force  coefficients  obtained  by  such  a  linearized  theory  are  generally 
correct  only  to  the  first  order.  In  the  present  rotational  development, 
terms  of  the  form  €y^  arise  in  the  boundary  conditions  for  the  per¬ 
turbation  velocities;  in  order  to  preserve  the  first  order  smallness  of 
these  velocities,  it  is  necessary  to  restrict  the  relative  vorticity 
€/u^  to  a  size  of  the  order  of  one.  In  this  case,  then,  the  flow 
reversal,  which  occurs  in  a  uniform  shear  flow  at  y  ***  U^€,  does  not 
occur  near  the  body-cavity  combination,  whose  ordinates  are  usually  very 

small  compared  to  one.  Finally,  from  the  Cauchy -Riemann  equations  as 

2 

applied  to  the  harmonic  perturbation  velocities,  it  is  seen  that  these 
velocities  change  very  slowly  in  space  when  the  streamline  slopes  and 
cuz^tures  are  small  [l] .  For  this  reason,  the  linearized  boundary 
conditions  may  be  applied  on  the  x-axls  instead  of  the  body-cavity 
surfaces.  These  boundary  conditions  for  the  harmonic  perturbation 
velocities  are  established  from  facts  known  about  the  body-cavity 
system  In  the  shear  flow. 

jja.  order  to  formulate  the  boundary  conditions  on  the  cavity,  it 
Is  first  necessary  to  define  a  cavitation  number.  The  cavitation 
number  is  the  parameter  which  relates  the  pressure  conditions  in  the 
fluid  stream  to  those  In  the  cavity.  It  Indicates  the  degree  of 
cavitation;  as  the  cavitation  number  decreases,  the  cavity  size 

p 

Increases.  Recalling  that  i|r  «  ^7  /2,  one  has  troa  equation 

(2.^),  to  the  first  order  in 


H  00  C 


-  1 


00 


00 


:s  constant. 


The  velocity  and  pressure  p^  are  taken  at  (-  oo,  o),  while  p^ 
and  q^  are  taken  on  the  cavity  streamline.  In  general,  the  extension 
of  the  cavity  streamline  does  not  coincide  with  the  x-axls  except  at 
the  leading  edge  of  the  slender  body.  The  stream  function  represents 
the  difference  between  the  harmonic  stream  function  value  on  the  x-axis 
and  its  value  on  the  up-stream  extension  of  the  cavity  streamline. 

g"’"  "  ■  ' 

See  equations  (2.13). 
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Thus,  ’IftT  *  •U  y  where  y  is  the  ordinate  of  the  cavity  streamline 

at  X  B  >ao.  In  the  case  of  a  uniform  flow  about  a  cylinder  with 

circulation  [l?],  -y  is  infinite.  However,  in  the  present  cases 

00 

of  a  linearized  uniform  shear  flow  past  slender  wedges  and  hydrofoils, 

the  value  of  y  is  not  known  a  priori.  It  follows  that  in  these 
00 

rotational  flows  the  cavitation  number  cr,  defined  in  the  usual  sense 
as 


u  = 


00  C 


(2.7) 


cannot  be  \ised  directly  since 


q  2ey 

^C  ,  oo 

a  =  -!r  -  1  + 


=2  2 

00 


00  00  00 

and  y  is  tmknown.  Therefore,  it  is  convenient  to  define  a  rotation- 

•'oo  ' 

al  cavitation  number  £.  Let  £  be  defined  so  that 


-  1; 


(2.7a) 


then,  in  the  following  treatment  most  terms  will  retain  the  same  form 
as  their  irrotational  counterparts.  When  either  €  -►  0  or  the  flow 
is  symmetric  £  e  a;  otherwise. 


2€y 


0  =  2  + 


OO 

00 


f 


This  relationship  provides  a  means  for  relating  the  pressure  difference 
and  flow  velocities  in  experimental  programs.  The  rotational  cavitation 
number  may  be  determined  either  by  measuring  the  up-stream  and  cavity 
velocities  directly  or  by  first  measuring  the  pressure  difference,  up¬ 
stream  velocity,  and  vorticity  and  then  measuring  y^,  the  ordinate 


■3 

For  a  symmetric  shear  flow  about  a  synmetrlc  slender  body,  the  x-axls 
and  the  bifurcated  cavity  streamline  do  coincide  for  x  0  and 
equation  (2.7)  holds.  In  this  case,  a  »  £.  See  Section  3.2  for 
discussion  of  such  synmetrlc  flow. 
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of  the  stagnation  streamline. 

On  the  cavity  vails,  from  equation  (2.^),  the  velocity  q  is 

c 

constant  because  is  constant  for  a  given  vortlclty  €.  Writing 

q  in  terms  of  the  perturbation  velocities  and  introducing  the  result 
0 

into  equation  (2.7a)  yields 


Z 


2u 

U 


2  2  2 

+  V  +  6  y^  .  2u€y^,)  . 


(2.7b) 


The  cavity  bovindary  condition  is  obtained  by  neglecting  terms  of  order 
greater  than  one  In  equation  (2.7b),  thus  producing 


U  Z 

=  u(x,0)  -  eyg(x),  a  <  X  <  i,  (2.8) 

which  is  applied  on  the  x-axls.  When  equation  (2.6)  is  applied  to 
the  upper  cavity  surface  in  the  flow  past  a  hydrofoil,  the  quantity  a 
equals  zero;  otherwise,  a  is  equal  to  one  for  flows  past  both  wedges 
and  hydrofoils.  If  one  lets 


U  =  U  +  u  -  €y 
2  p  2 

on  the  cavity,  then  q  *=  U~  +  v  .  Thus,  from  equation  (2.7b),  to  the 

0  0 

first  order, 

(1  +  2/2).  (2.9) 

Fr<xn  equation  (2.7a)  one  obtains 

U  «  q  B  U  Vl  +  2  ^ 

C  C 

In  the  present  work  equation  (2.9),  which  is  consistent  with  the 
linearisation,  is  used;  however,  it  should  be  noted  that 

1  +  z/2  «  Vl  +  2  '. 


■*1 
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It  was  noted  previously  that  the  harmonic  perturbation  velocities 
change  very  slowly  when  the  streamlines  have  small  slope  and  curvature. 
For  this  reason^  since  the  cavity  Is  assumed  to  be  long  and  slender,  the 
variations  in  the  quantity  ey  in  equation  (2.8)  are  certainly  small 

over  most  of  the  cavity.  It  Is  reasonable,  then,  to  replace  €y  with 

+  -  ° 
an  average  value  -  e  on  the  upper  and  lower  cavity  surfaces  respective¬ 
ly.  This  averaging  technique  was  Introduced  by  Parkin  [U]  In  dealing 
with  a  linearized  theory  for  flow  past  a  hydrofoil  In  a  transverse 
gravity  field.  The  averaged  quantity  Is  non-dimensional.  It 

must  be  treated  as  part  of  the  problem's  solution  and  Is  used  as  a 
basic  parameter  when  the  resxilts  cannot  be  determined  directly  as  a 
function  of  ®/u  .  Note  that  the  relative  vorticlty  e/u  ,  which  has 
the  dimensions  (l/lENQTE),  Is  Independent  of  the  length  of  the  solid 
body  and  truly  measiares  the  relative  size  of  the  small  constant  vor- 
tlcity  and  the  up-stream  velocity,  l.e.,  the  rate  of  change  of  velocity 
with  y  at  infinity.  Equation  (2.8)  now  becomes 

u(x,0)  *=  U^2/2  +e, 

and  u(x,0)  =  U^E/2  -  i,  1  <  x  <  i,  y  <  O.J  (2.8a) 

These  equations  are  applied  on  the  upper  and  lower  cavity  surfaces 
resi>ectively. 

The  total  velocity  on  the  solid  body  must  be  tangent  to  the 
surface  of  the  bodyj  hence,  in  terms  of  the  perturbation  velocities, 

v(x,y^) 

^  ’ 

The  denominator  may  be  written  as 

with  Tj  m  (u(x,y  )  -  u(x,y  )  -  «y^(x)  +  6y_(x)).  From  the  basic 

a8sunq)tlons  of  the  theory,  t)  <  <  u  •  Since 

c 


I 
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V 


_ 

the  denominator  may  be  e:q>anded  in  a  binomial  series,  ^e  res\ilt  is 


V  r, 


After  linearization,  the  bovindary  condition  is 


v(x,0) 

— ST"  =  “tP^  » 

C 


0  <  X  <  1  , 


(2.10) 


and  it  is  applied  on  the  x<axls.  Furthermore,  this  equation  is  also 
valid  on  the  cavity  surfaces  and  gives  the  surface  slope  at  any  point 
along  the  x-axls.  In  the  cases  of  slender  wedges  or  hydrofoils, 

clyJx) 

=  teui  a  te  a  . 


Equation  (2.10)  becomes 

v(x,0)  =  +aU^,  0<x<l,  (2.10a) 

on  the  upper  and  lower  wedge  surfaces  respectively,  and 


v(x,0)  »-aU,  0<x<l,  (2.10a) 

C 

on  the  solid  surface  of  the  hydrofoil. 
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Examination  of  equations  (2.8a)  and  (2.10)  shows  that  the  result 
acquired  in  equation  (2.10)  provides  automatically  for  smooth  separation. 
This  condition  may  be  thou^t  of  as  equivalent  to  the  Kutta  condition  in 
airfoil  theory  in  that  both  conditions  serve  to  single  out  a  vmlque 
solution  to  the  flow  boundary  value  problem. 

The  condition  of  cavity  closure  is  characteristic  of  Tulin's  model 
of  the  finite  cavity.  If  the  rotational  and  new  irrotational  flows 
are  studied,  it  is  seen  that  the  subtracted  portion  \|fp  of  the 
rotational  flow  makes  no  net  contribution  to  the  flow  within  the 
cavity-body  shape.  Hence,  the  cavity  closure  condition  must  hold  in 
either  flow,  l.e.,  the  net  strength  of  sources  within  the  body-cavity 
system  must  be  zero. 

Finally,  since  the  base  flow  is  tindlsturbed  at  infinity,  the  per¬ 
turbation  velocities  (u,v)  approach  zero  at  great  distances  from  the 
body-cavity  system.  In  sumnary,  the  conditions  to  be  itqposed  are: 


U  L 

u  ■  ®  ^be  cavity 

V  =  +  0!  U  on  the  body 
—  c 

(u,v)  -►0  at  infinity 
The  cavity  is  closed. 

The  flow  must  separate  smoothly  from  the  solid  body. 


These  conditions  are  sufficient  to  determine  the  harmonic  perturbation 
velocities  (u,v). 

2.2.  Methods  of  Solution 

Two  methods  of  solution  and  their  corresponding  boundary  value 
problems  are  given.  The  first  is  that  method  used  originally  by  Tulin  [l]. 
This  method  makes  use  of  the  velocity  potential  of  a  sovirce  distribution 
and  singular  integral  equation  theory.  The  second  method  is  based  on 
conformal  mapping  of  the  physlceJ.  plane  onto  the  exterior  of  the  unit 
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circle.  This  latter  method,  which  has  been  used  by  .others  [4,  Q,  l8], 
has  proven  most  convenient  for  the  present  work;  however,  the  former 
method  Is  outlined  here  for  completeness  and  comparison. 

The  first  method  Is  given  here  for  a  symmetric  flow  such  as  that 
found  In  Section  3*2  and  Figure  4.  The  solution  for  the  mathematical 
problem  arising  from  the  given  conditions  on  (u,v)  Is  found  In  terms  of 
the  velocity  potential  which  exists  for  the  harmonic  flow  represented 
by  (u,v)  and  The  potential  of  a  distribution  of  sources  along 

the  x-axls  between  the  leading  edge  of  the  body  and  the  close  of  the 
cavity  Is 

i 

0(x.,y)  «  -  in  r  d|  .  (2.11) 

The  distribution  function  4(1)  represents  the  source  strength  and 

2  /  \2  2 
r  -  (x  -  1)  +  y  . 


This  distribution  produces  a  function  with  the  required  symmetry. 
The  velocities  are  designated  by 


u(x,y) 


-^4 


aii) 


i2Lzjl 


ki(5)  ^ 


(2.12) 


The  given  boundary  conditions  on  (u,v)  are  applied  on  the  x-axls,  and 
one  must  pass  to  the  limit  y  0  In  these  equations.  Taking  the 
Oauchy  principal  value  of  the  Integrals  when  required  and  using  the 
substitution  tan  t  ■  (x  -  |)/y,  one  obtains 


u(x,0) 


1 

31 


v(x,0)  «  n  (x)  /  2 


(2.12a) 


Incorporating  equations  (2.12a)  Into  the  conditions  given  at  the  end 
of  Section  2.1  produces  the  following  complete  boundary  value  problem: 
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To  find  ^(x)  such  that 


^  1  dy^Ct) 

/  .  mjj:  +  an?  -  f  1?^" '  dt  .  f  (x),  i  <  x  <  i. 


ciy«(3c) 

b.  n(x)  -  2U^  -g2 - ,  0  <  X  <  1. 


c.  The  cavity  Is  closed;  hence 

I 

fi(x)dx 


dx  -  -  2U  y  (0). 


C'O 


d.  The  separation  Is  smooth;  hence  ^(x)  finite  as 

X  ^  X  a 


Since  the  flow  Is  synmetrlc,  it  Is  sufficient  to  consider  only  y  >  0. 
From  equation  (A-^)  of  Appendix  A,  one  has 


Ia(x) 


X 

n^Y(i  -  x)  (x 


Once  ^(x)  Is  known,  the  cavity  shape  and  pressure  force  coefficients 
can  be  calctilated  by  using  equations  (2.12a)  and  the  definitions 
given  m  Section  2.3* 

The  second  method  of  solution  also  makes  use  of  the  new  harmonic 
flow.  Ih  this  case.  It  is  convenient  to  work  with  ccasplex  variables. 
FoUovlng  substitution  of  (u,v)  and  into  the  vortlclty  equation 
(2.1)  and  continuity  equation  (2.2)  they  reduce  to 


and 


bv  &u " 
du  -dv 


(2.13) 


respectively.  In  terms  of  the  complex  variable  z  ■  x  +  ly,  the  total 
ccssplex  velocity  may  be  defined  as 


W(z)  »  +  i«(2  -  z)  +  w(z). 


I 
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with  w(z),  the  complex  perturhatlon  velocity,  given  by  w(z)  ■  u  -  iv. 
Equations  (2.13)  are  seen  to  be  the  Cauchy-Riemann  equations  for  w(z); 
from  these  equations  and  the  continuity  of  the  flow,  it  follows  that 
w(z)  is  analytic  outside  the  cavity-body  system.  Also,  since  (u,v)  -*■  0 
at  infinity,  w(z)  must  vanish  at  infinity. 

The  boundazy  conditions  are  applied  along  a  slit  in  the  z-plane 
corresponding  to  the  x-axis  where  0  <  z  <  1  and  z  is  real.  The 
complex  velocity  w(z)  is  analytic  in  the  region  exterior  to  the  slit 
in  the  z-plane.  The  complex  z-plane  and  typical  boundaxy  conditions 
are  shown  in  Figure  3>  The  mathematical  problem  is  to  ascertain  an 
analytic  function  w(z)  subject  to  boundary  conditions  on  both  its  real 
and  imaginary  parts,  as  given  at  the  end  of  Section  2.1. 

The  specific  problems  for  the  wedge  and  hydrofoil  are  tabulated 
in  Table  1.  Note  that  in  this  table,  conditions  e.  and  f.  are 
equivalent  to  requiring  that 

Ln(w)  B  0(i),  z  -►  « 

Real(w)  B 

The  conditions  listed  in  Table  1  are  sufficient  to  determine  w(z).  It 
is  expected  that  the  complex  perturbation  velocity  w(z)  will  exhibit 
singular  behavior  at  those  points  on  the  slit  where  large  changes 
occur  in  the  magnitude  and  direction  of  the  velocity.  These  large 
changes  occur  at  the  leading  edge  of  the  solid  body  and  at  the  closure 
of  the  cavity.  Condition  g.  in  Table  1  is  imposed  to  restrict 
the  magnitude  of  the  singular  behavior  so  that  the  pressure  distribu¬ 
tion,  which  is  proportional  to  Re(w),  remains  Integrable. 

The  determination  of  w(z)  is  accomplished  by  a  method  of  conformal 
mapping  which  permits  the  use  of  resvilts  edready  known  in  airfoil 
theory.  In  accordance  with  Wu  [l8]  and  Parkin  [4],  the  z-plane 
is  mapped  conformally  by  a  succession  of  transformations  onto  the 
C-plane.  The  complex  velocity  w  is  held  invariant  at  corresponding 
points  of  the  mappings.  The  transformations  are  listed  in  Table  2, 
and  the  varloxis  mapped  planes  are  shown  in  Figures  3  and  in  which 
corresponding  regions  and  boundary  values  are  shown. 
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It  l8  seen  that  the  houndary  value  problem  for  v  may  be  solved 
In  any  of  the  transformed  planes.  The  C-plane  is  used  directly 
for  ved^^s,  while  the  problem  is  solved  in  the  Q-plane  for  hydrofoils. 
In  either  case,  however,  the  singular  functions  which  are  the  basis 
of  the  problem's  solution  are  established  in  the  C** Plane.  From 
Figures  3  and  5,  it  is  clear  that  the  real  parts  of  the  singular 
functions  must  satisfy  particular  conditions  on  the  real  axis,  while 
the  Imaginary  parts  must  satisfy  other  conditions  on  the  unit  circle. 

In  addition,  the  limiting  conditions  on  w  as  z  «  must  be 
met  at  the  corresponding  points  in  the  C'Plana  or  Q- plane.  A 
complete  solution  w(^)  is  formed  from  a  series  of  singular  complex 
functions.  These  functions,  familiar  in  airfoil  theory,  and  their 
proi>ertles  are  listed  in  Table  3>  The  strength  of  the  singularities  is 
limited  by  the  pressure  integrabllity  condition  previously  noted. 

The  singularities  are  not  of  higher  order  than  simple  poles  at  points 
on  the  slit.  Figure  3  shows  that  the  leading  edge  of  the  hydrofoil 
maps  onto  C  1;  the  leading  edge  of  the  wedge  onto  C  ■  cuid  the 
cavity  closure  onto  C  *  These  are  the  centers  of  the  singular 
behavior  in  the  linearized  theory;  hence,  most  of  the  v^(Z)  in 
Table  3  are  singular  at  these  points.  These  functions  in  Table  3  have 
already  been  extensively  used  in  linearized  cavity  theory  and  are 
discussed  in  some  detail  by  Parkin  [4,  6]  and  WU  [16,  19] . 

The  solution  w({|)  to  a  particular  problem  is  given  in  terms  of 
the  singular  functions  w^(C)  by 

w(C)  +  M  +  IH, 

1 

where  the  constants  K^,  M,  and  N  are  eissumed  to  be  resQ..  These 
constants  are  then  determined  by  the  conditions  given  in  Table  1, 
since  the  w(C)  formed  is  already  analytic  off  the  slit.  In  addition 
to  establishing  the  constants,  the  given  conditions  s^jso  pzx)duce  a 
relationship  between  the  cavitation  number  £  and  the  cavity  length  i, 
with  €  and  a  as  parameters.  Following  determination  of  w,  the 
cavity  shape  auid  pressure  force  coefficients  can  be  found  by  using 
equation  (2.10)  and  the  definitions  given  in  the  following  paz*t  of 
this  section. 
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2.3*  Calculation  of  Resiilts 

The  results  of  the  linearized  theory  include  the  length- cavitation 
numiber  ratio  (already  found  In  Section  2.2),  the  cavity  shape,  and  the 
pressure  force  coefficients  for  lift  (C^^),  drag  (Cjj),  and  moment 
(Cjjq).  From  equation  (2.10),  the  cavity  shape  is  given  by 

®  (2.Ht) 

For  wedges,  a  =  1  on  both  surfaces,  while  yQ(l)  ■  +  a  on  the  upper 
and  lower  surfaces  respectively.  For  hydrofoils,  a  «  0  and  7^(0)  ■  ( 


on  the  upper  cavity  surface;  a  =  1  with  yQ(l) 


a  on  the  lower 


surface.  Also,  using  equation  (2.10),  one  may  calculate  the  body- 
cavity  area  S,  which  Is 


This  equation  may  be  Integrated  by  parts  and,  since  the  cavity  must 
close,  l.e.,^dy  ^  0,  one  has 


vxdx. 


(2.15) 


The  contour  Integral  B-K!  follows  a  closed,  counter-clockwise  path 
over  the  surfaces  of  the  body-cavity  combination.  Recalling  that 
V  se  -Im  w(z),  one  may  write  the  previous  equations  so  that 


ycW 


wclz  +  yo(a) 


(2.lUa) 


'k^J 


wzdz. 


(2.15a) 


where  x  s  z  on  the  slit  In  the  complex  z-plane.  Tlus  latter  result 
was  first  given  by  Geurst  and  Tlinman  [20]. 

Calculation  of  the  pressure  force  coefficients  Is  based  on  the 
pressure  coefficient  C^,  which  is  given  by 
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(2.16) 


In  each  flow  the  solid  hody  and  cavity  surfaces  lie  on  the  sane  stream¬ 
line.  Thus,  the  BemoxiUl  equation  (2.^)  yields 

P  -  (1  -  %/%)• 

The  term  on  the  right  In  this  equation  nay  he  simplified  by  lineariza¬ 
tion.  Writing  q^  ■  and  qb  «  +  v^,  one  has  q^  =  to  the 

c  C  O  '  c  c 

first  order  and 

1  -  -  -^(qo  -  -  ('lo  • 

-  -2(0,,  -  U^)/U,  -  (u^  - 

-  -2(0„  -  Je)/»c 

after  linearization.  The  pressure  coefficient  may  nov  be  written  as 

®p  -  -20c<"o  - 

with  and  being  the  x-conq>onents  of  velocity  on  the  solid 
body  and  cavity  surfaces  respectively.  From  eqxiatlon  (2.6a), 

Uq  =  u(x,y^)  ”  «3^  +  U^,  and  from  equation  (2.9) 

Thus, 

Cp  -  -  2(1  +  |)  (u(x,y^)  -  6y^-  Uj)/U^.  (2.16a) 

There  has  been  considerable  discussion  [U,  6,  19  ]  x^g^trdlng  the 

appropriate  method  of  defining  and  linearizing  the  pressure  coefficient 

so  as  not  to  lose  any  first  order  terms.  The  above  result  Is  consistent 

with  that  given  by  Wu  D-9  1/  in  that  the  flow  Is  continuous  and  C  s  0 

P 

at  the  trailing  edges  of  the  solid  body  for  €  ■>  0.  As  Parkin  has 
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pointed  out  [U],  however,  in  order  to  he  consistent  with  the  averaging 
approximation  when  €  /  0  and  to  have  a  continuous  velocity  at  the 
trailing  edge,  one  must  take  [eyl  b  |e|  on  the  solid  body  for  x  *  1. 
On  a  wedge,  y  =  +  ox;  so,  to  the  order  of  the  approximation, 

ey  =  +  €  X 

on  the  upper  and  lower  wedge  surfaces  respectively.  On  a  hydrofoil, 
y  B  -  OK,  so 


€y  B  -  6x 


on  the  foil  to  the  order  of  the  averaging  approximation.  The  final 
result  for  the  pressure  coefficient  is  given  in  Table  4.  In  addition, 
this  table  gives  the  remaining  coefficients  in  terms  of  Cp.  Note 
that  the  drag  coefficient  in  the  wedge  flows  is  based  on  the  base 
area  of  the  wedge;  the  remaining  coefficients  use  the  unit  chord  of 
the  body  as  a  characteristic  length.  The  method  used  to  calculate 
these  force  coefficients  follows  that  used  by  Wu  [l8]  on  wedges 
and  Parkin  [4]  on  the  hydrofoil. 

Finally,  a  rational  means  for  evaluating  the  average  parameter 
e/u^  must  be  determined.  From  the  manner  in  \rtilch  5/u^  arises, 
it  is  mathematically  reasonable  to  let 

^  (cavity  area)/2(i  -  l), 

l.e.,  e/u^  equals  the  average  value  of  ^ly^,!  over  the  cavity.  On 
the  other  hand,  the  vorticlty  c  creates  an  additional  circulation 
in  the  flow,  and  it  is  this  important  flow  property  which  characterizes 
the  influence  of  the  vorticlty.  (This  result  holds  also  in  the  case 
of  a  supercavltatlng  flow  past  a  hydrofoil  in  a  transverse  gravity 
field  [4]).  Thus,  is  chosen  so  that  the  actual  clrcxilatlon  F 

is  equal  to  the  circulation  F  based  on  the  constant  perturbation 
velocities  associated  with 
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The  circulation  P  Is  defined  as 


da 


where  q  Is  the  total  vector  velocity  and  ^  Is  the  elemental 

vector  path  length  on  the  closed  contotir  where  P  Is  measured.  Since 

the  perturbation  velocities  are  defined  In  the  harmonic  rather  than 

the  physical  z-plane.  It  Is  convenient  to  balance  P  and  P  In 

the  harmonic  plane.  From  equations  (2.6)  and  (2.6a),  one  finds  that 

U  »  U  +  u  -  6y  Implies  that  U_  «  U  +  u  and  V  ■  v  inplles  that 
00  £1  00 

V„  «  V.  Therefore,  on  the  cavity  where  U  *»  U  , 

XI  c 


+  ey 


and  on  the  solid  bodies. 


Vg  -  T  . 


The  actual  circulation  P  and  T  are  given  to  the  first  order  In 
Table  Note  that,  to  the  first  order  on  the  solid  bodies,  the 
circulation  integral 


+  u)  dx  +  vdy] 


+  u)dx. 


body 


body 


Observe  also  that  In  the  e3q)re8slon  for  in  Table  the  terms  in 

the  large  parenthesis  are  precisely  equal  to  the  cavity  area  for  the 
wed{^  flow  and  equal  to  the  cavity  area  less  the  triangular  area 
between  the  hydrofoil  and  the  x-axls  for  the  hydrofoil  flow.  Thus, 
the  matching  of  clrc\ilatlons  at  once  provides  a  rigorous  and  Intuitively 
satisfying  result. 
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3.  APPLICAIIOIE  OF  THE  LIMARIZED  THEORY 


In  this  section  the  linearized  theory  developed  In  Section  2  Is 
applied  to  three  problems :  the  first  two  are  the  supercavitatlng, 
asjonmetrlc  and  symmetric  shear  flows  past  slender  wedges;  the  third 
concerns  a  supercavltatlng,  uniform  shear  flow  past  a  flat-plate 
hydrofoil.  The  solutions  to  all  three  problems  are  found  by  means  of 
conformal  mapping,  using  the  singularities  shown  In  Table  3. 

3.1.  Asymmetric  Flow  Past  a  Wedge 

This  first  flow  Is  a  supercavltatlng,  uniform  shear  flow  past 
a  wedge  of  unit  length.  The  asymmetric,  undisturbed  velocity 
profile  of  the  flow  Is  shown  In  Figure  1  (solid  line  profile).  The 
notation  used  Is  that  given  In  Section  2  and  Figure  2a. 


3.1.1.  Solution  of  the  boundary  value  problem 

The  mapping  of  the  physical  z-plane  and  the  appropriate  boundary 
cozxdltlons  are  Indicated  In  Figure  3&*  From  the  conformal  trans¬ 
formations  In  Table  2,  It  Is  easily  shown  [^]  that 


[  («'  -  C^)  J 

where  and  C2  roots  of 

C**  +  ac®(2i  -  1)  -  (C^  -  C^)  (C®  -  ^) 


0. 


These  roots  are 


■  1  ( •/T  +  V/  -  i)  ^ 
Cg  - 1  (VT-  vrri). 


(3.1) 


(3.2) 


It  Is  seen  from  Figure  3a  that  Is  outside  the  unit  circle  and 
represents  the  point  z  ■  m,  while  Is  Inside  the  unit  circle  and, 
hence,  does  not  represent  a  point  of  the  physical  plane. 
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The  cooplete  houndary  value  problem  for  this  wedge  flow  Is  given 
In  Table  1,  which  was  daveloped  In  Section  2.  By  using  the  singular 
functions  listed  In  Table  3,  It  Is  possible  to  constiruct  a  solution 
function  w(C).  A  comparison  of  the  boundary  conditions  and  available 
singularities  shows  that  w(C)  should  have  the  form 

w(0 - ^  in  ^  +  1A(C  -  +  B 

+  ICinC  +  ID  ^  .  (3.3) 

The  function  l(C  -  l/C)  is  selected  to  provide  the  proper  closure 
singularity,  and  the  tenu  l(^  •  i)/(C  i)  is  used  to  satisfy 

the  condition  that  w(z)  0  as  z  The  remaining  functions  are 

needed  to  fulfill  the  boundary  conditions  (a)  through  (d)  in  Table  1. 
The  real  constants  B  and  C  are  determined  by  considering  w(C)  on 
the  real  axis  of  the  ^‘plane*  Equating  the  values  of  w(C)  and  the 
bomdary  conditions  there,  one  has 

U  Z 

w(C)  “  B  «  +  ®,  >  0 

- 

w(5)  »  B  •  CII  »  g '  •  C  ^ 

Therefore,  B  ■  U^2/2  +  e  and  C  **  2«/lI.  Equation  (3.3)  can  now 
be  rewritten  as 

U  Z  2DtU  *  . 

'^(C)  “  +  €  -  ||.  ^  in  +  i 

+  1A(C  -  i)  +  ID  Vl-i  .  (3.3a) 

^  C  +  1 

In  this  form  w(C)  satisfies  the  differential  equation  and  conditions 
(a),  (b),  (c),  (d),  (g),  and  (h)  of  the  boundary  value  problem.  The 
remaining  conditions  will  serve  to  establish  A  and  D  axid  to  determine 
a  unique  relation  between  i  and  Z  for  given  values  of  s  and  a. 
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3*1*2.  Results 

As  Wu  has  pointed  out  [l8],  much  can  be  learned  by  expanding  v(z) 
in  negative  powers  of  z  as  z  ->  ««  The  result  Is  an  Infinite  series 
of  the  form 


w(z)  «  a  +  lb  + 
o  o 


a^^  +  Ibj^  ag  +  Ibg 


+  0 


(3.4) 


From  botindary  conditions  (e)  and  (f),  one  has 


*  a^  “ 


0. 


As  z  ->  00,  5  hence,  using  equations  (3*1)  and  (3.2),  C  can 

be  expanded  In  descending  powers  of  z.  Acosta  [^]  has  shown  that 
one  obtains 


c 


«1 


li  +  *  vg(7TTr  t2i  - ./itJTT)] 


8z‘ 


(3.5) 

This  result  Is  now  Introduced  Into  equation  (3.3a)  and  the  coaft>lDatlon 
Is  simplified.  This  process  Is  accoiqpllshed  In  several  steps.  First, 
from  Acosta's  work: 


U  L  _  2c«J^  .  ^  T  /  U  Z  M  , 

-5-  +  e  -  -j-  in  +  U(C  -  f)  -  (-5-  +  «  -  -5-  in  ^  ^ 

.2Avr)+i  [-jS  VT.  A(i - 1) /r] + ^  (i  +  i)-/r 

z 

“  ^  (“^  "  1)  (3-®  +  1)  • 

Second, 

l^cin  (|b1^c  [inC^^  +  Ai  (1  Q)], 

with  Q  -  h/z  +  h/z^  +  0(l/z^),  M  «  Vi(i  -  1)’  /2,  and  H  -  Vi(i  -  1)* 
(21  +  1  +  •JKjl  -  l)*)/8.  Following  a  series  expansion  on  the  small 
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parameter  Q,  one  has 


K  -  -»  *  #  H-/T*  ^/^3)  *  *  0^^) 


Finally, 


r2  .  1  Cl  (1  +  Q)^  -  1 

ID  K - i  -  ID  -i - n - , 

r  +  1  CT(1  +  +  1 


with  Q  as  hefore.  Letting  Ci  ®  “  (2i  -  1  +  2Vi(i  -  l))  *  -  T,  one 
obtains 

IP  ^  ir  ^  ^  , 

C^  +  1  (1  -  T)  -  T(2Q  +  Qi) 

and  by  setting 


b  «  (1  +  T)/T,  a  -  (1  -  T)/T,  and  n  «  2Q  +  QT, 


he  has 


for  T)  «(a,b).  By  expanding  the  denominator  of  this  result, 
neglecting  terms  of  order  (l/z^),  and  simplifying  the  result,  one 
obtains  finally 


iD 


rTi 


T  .  UT  M  .  UMT  1 

^  (1  -  T)‘  ■ '  (1  -  if  ■  7 


The  series  e:q>ansion  of  v(b)  in  descending  powers  of  z  is  a 
combination  of  these  resxilts,  and 
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W(z)  =  ^  in  +  l|i  in(>/]r+  ^/^^l)  -  iD^-f|j 


z 


(1-T)‘  j 


z 


ou 


(ifi)  sir  -  ^  (i-i)  (3^1)  v/^+  In  ®  (2^1) 


-  ID- 


+  0 


(1-T)- 


C?) 


+  s/i(i-l))(l-T)  +  i(i-l)  ^ 


(3.6) 


Since  It  Is  required  that  «  h^  b  0,  It  follows  from  equation  (3.6) 
that 


A  «  0U^/ll(i-l) 


(3.7a) 


and 


®  ®  1^  M  >/^+  s/i  -  1). 


(3.7b) 


Using  these  results  and  the  further  requlretoent  that  a^  >  0,  one  has 


U  2  0«J^ 

«o  C 

“2 — r 


Because  U  ■  U  (l  l/2),  the  above  relation  can  he  simplified 
c  ®® 


to  give 


r  a/«_*/r+i  2Jr\ 


(3.8) 
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This  equation  relates  the  cavitation  nuniber  L,  the  cavity  length 
and  the  wedge  semi-angle  a.  It  is  seen  to  he  the  same  result  as  that 
obtained  hy  Tulin  [1]  and  Va  [l8]  for  an  Irrotational  wedge  flow. 

In  this  case  the  rotationallty  of  the  flow  has  no  first  order  effect 
on  the  cavity  length.  This  same  conclusion  can  he  reached  hy  purely 
physical  reasoning  as  follows:  if  a  first  order  length  effect  enters 
as  Ks,  then  changing  the  sign  of  a  reverses  the  first  order  effect 
hut  simply  Inverts  the  flow  field.  Hence,  a  contradiction  vo\;ild  result 
from  the  presence  of  a  first  order  length  effect.  This  same  reasoning 
applies  to  first-order  changes  in  the  drag  and  cavity  area,  hut  not  to 
changes  in  the  other  pressure  force  coefficients  or  the  cavity  shape. 

As  a  result  of  the  above  information,  the  complex  perturhatlon 
velocity  v(z)  may  now  he  written  as 


i^re  it  is  seen  that  a  is  the  area  of  the  wedge.  Because  the 
complex  velocity  w(z)  is  azutlytlc  off  the  slit. 
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B+C  T 

where  T  Is  a  circle  of  large  radius  (|z|  ->00)  surrovindlng  the  slit. 
The  slit  and  contour  path  lines  are  plotted  In  Figure  6.  Applying 
the  theory  of  residues,  one  has 


and  from  equation  (3 >9) 


•*  SHag, 


*^2 


=  -  oaj^i2/^/2n. 


Th\is, 


S  =  (3.11a) 


To  find  the  cavity  shape,  one  rewrites  equation  (2.lUa)  as 

yc-a-lm^  2i^gdC 
on  the  iqpper  cavity  surface,  and 

on  the  lower  surface.  On  the  cavity  surfaces  (see  Figure  3a),  both  C 
and  z  are  real  and  x  ■  z.  The  cavity  ordinate  y^  Is  related  to 
the  cavity  abscissa  x  by  determining  x  as  a  function  of  the  para*> 
meter  t.  The  derivative  dz/dC  Is  found  by  dlffeirentlating  equation 
(3.1),  and 


(3.12) 


dz  8i(i  -  i)c(r  - 1) 


^  (r  +  T)^  (C^  +  R)^’ 


The  relations  Cg  -  1>  »  -T,  aiid  Cg  ■  "R  -  “t2i  -  1  -  2 

have  heen  used  to  achieve  this  result.  When  C  Is  real,  the  imaginary 

part  of  v(^)  in  equation  (3.3&)  becomes 


l40«„ 

Im[w(C)]  =  - 


tan“^  ^  ^  «  AilCl  +  A(C  -  -J)  +  D  ‘ 


Combining  the  above  results  and  vising  equations  (3.7)^  one  has  on 
the  cavity  surfaces  the  following: 


a.  on  the  upper  surface,  t  >  1, 


a  -  jijut  <c-  t>  *  5  (i^  ■  ^ 


b.  on  the  lower  surface,  t  <  -1, 


(3.13a) 


y  ■  -a  +  8i(j8-l) 


a  /,  1,  IS  l-X  info/T+.^n) 

flTyK  -  v)  -  BOTTS  •  '  lili)  '  ^ 


4c 

’  Btr  Si 


•1 1 1 .  ittkia!;, — ,• 


(3.13b) 


c.  from  equation  (3.1)#  for  t  >  1, 


ti" . 

L  (t2-KP)(t2+R)_ 


(3.13c) 


Next,  the  pressure  force  coefficients  are  determined.  The 

pressure  coefficient  C  was  developed  in  Section  2.3  and  is  given 

P 


I 
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explicitly  In  Table  4;  the  other  coefficients  are  given  In  terms  of 

C  .  From  Table 
P 

Cp  =  (2+E)  r  +  +  I  -  g-l  .  (3.14) 

U  00  ^  aJ 

On  the  wedge  surfaces^  C  e  .  The  upper  surface  corresponds  to  the 
plus  sign  In  the  Cp  expression  and  0  <  0  <  II/2,  while  the  lower 
surface  corresponds  to  the  minus  sign  and  II/2  <  d  <  H.  On  the  wedge^ 
0  <  X  <  1.  The  combination  of  equations  (3.3a)  and  (3.7)  and  the 
Introduction  of  ?  »  e  gives  the  complex  velocity  on  the  wedge 


v(e)  E  /o  X  “ 

U  »  -5  -  (2  +  2)  ^ 


in 


e^®  +  1  1 


/  Id 

(e  -  e 


■^)1, 


120 


^  r  m^^+a/TTi)  ‘  +  5"  (1 "  Tt0)- 

00  e  +  1  » 

Simplifying  this  resvilt  and  taking  the  real  pairt  produces 


sin  0  -  Y  7'  Aa(./r  +  Ji  •  1)  tan  0 

^00 


+  T  ^  "  (2  +  E)  §  /n  rv  . 

00  '  I 

From  the  wedge  flow  transformations  listed  In  Table  2, 

i  ({ +  -  (i  - 1)  . 

16 

Since  x  «  Re  (z)  and  d  s  e  on  the  wedge,  one  obtains,  after 
scxne  simplification, 

^  /  cos^0  .  (3.15) 

i  -  stn^0 

The  above  resxilts  are  now  introduced  Into  equation  (3*14)  and 
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(2  +  2:)|§  (2  +  z)  (^§-f  +  in  |i 

#• 


1  +  Bln  e 
sin  d 


I) 


6 

•^u 


2d 

TT 


-  1  +  wfeil  tan  e  .  in  (/T  +  /Tl)  + 


i  cos^  e 

i-sln  6 


J 

(3.16) 


The  positive  sign  and  0  <  d  <  Il/2  are  used  on  the  upper  wedge  surface 
and  the  negative  sign  and  II/2  <  d  <  n  on  the  lower  sturface.  Equations 
(3.13)  and  (3.16)  then  give  the  pressure  coefficient  as  a  function 
of  X  on  the  wedge  surfaces. 

The  drag  coefficient  based  on  the  baise  width  of  the  wedge  is,  from 
Table  4, 


C 


D 


(3.17) 


The  contour  Integral  V  follows  a  closed,  counter-clockwise  path  on 
the  wedge  surfaces.  If  equation  (3.1^)  Is  Introduced  into  (3.17)^  and 


the  relationships 


dy.— 4x 

c 


1 

55 


jn 


$5 


and  2uv 


dz  +  (2+E) 


-  Mv®) 


are  used,  then 


^(7^x.f)ay 


, (3.17a) 


since  dx  *  dz  on  the  slit.  The  first  quantity  in  the  brackets  above 
can  be  written  in  terms  of  the  contour  Integration  paths  shown  in 
Figure  6  because  w  is  analytic  off  the  slit.  Hence, 


^w^dz  -^w^dz  -  ^ 


>w^dz. 


The  path  T  Is  a  large  radius  circle,  as  before;  the  path  C  around 
the  cavity  on  the  slit  consists  of  the  cavity  walls  pl\is  a  small  circle  e 
(radius  t  •*  O)  which  sxurounds  the  point  z  ■  i.  On  T,  w(z)  Is  of 
the  form 
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so  that 


ih^  +  Ih- 


+  0(h), 


w^(z)  =  -  -|  +  0(ir). 

Z  Z'^ 


Thus,  hy  the  theory  of  residues. 


^w^dz  B  0. 

T 

The  Integral  over  the  cavity  Is  given  hy 


2. 

V  dz 


C 


a^4y 


•/(%*nj./2)„U^dy 


4-  £n  J, 


T# 


s  J 

2 

where  Im  v  «  -  2uv  on  the  silt,  and  Is  the  Integral  over  e. 

The  Integrals  S  over  the  silt  are  taken  on  the  lower  and  upper  cavity 
surfaces  respectively  as  Indicated  hy  the  stjhscrlpt  on  the  Integrand. 

On  the  cavity  €y  «  +  €,  on  the  wedge  +  "ix  *  ey  «  ®y/c<.  It 

C  o 

follows  from  the  definition  of  a  closed  contour  IntegreO.  that 


^^^ydy  ■  0  and  €dy  *  0. 
W  S 

Thus,  equation  (3.17&)  in^y  he  simplified  to 


(3.17b) 


In  order  to  evaluate  this  resxilt,  w(z)  must  he  eiqpanded  as  z  •*  i. 
Taking  equation  (3.1)  and  letting  z  ■*  t  gives 
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g  ^  as  i. 

JTTl 

The  expansion  of  v(z)  about  z  m  £  and  evailuatlon  of  the  resulting 
Integrals  over  the  e  path  are  accomplished  In  Appendix  B.  There  It 


Is  found  that 


-^w^dz  «  8II  -  1). 


Introducing  this  result  and  the  value  for  A  given  In  equation  (3«7&) 
Into  eqmtlon  (3.17b)  above  gives  the  drag  coefficient 


a(2  +  Lri 

gfl  .  iV  • 


(3.18) 


As  expected,  this  result  is  precisely  that  given  by  Tulin  [l)  and  Wu  [18] 
for  the  Irrotatlonal  flow  aXtout  a  wedge,  end  the  rotation  has  no  first 
order  effect  on  the  drag. 

By  utilizing  the  previously  detexvlned  results  of  Table  U,  the 
lift  coefficient  is  found  to  be 

-  -(2  +  dx  +  (2  +  E)^  ^  x)  to. 

W  **  W  * 

The  real  part  of  v(z)dz  Is  equal  to  uto  on  the  wedge  and 


-  -  (2  +  E)  Re  ^5-  to  -  (2  +  E)  ^ 
From  equation  (3*9)  theory  of  residues 


wto  -  2ni(ibj^), 


(3.19) 
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since  a.^  »  0,  and  jfo(^/r+/s/E^)j  . 

Recalling  that  w  is  analytic  off  the  slit,  one  may  write 


^wdz  =^wdz  -^wdz. 
W  T  C 


Then, 


Re  ^wdz 
W 


and 


^wdz  =  -2€ViK(i-l)'  ^1  +  in  +  Vi-l)'J  -  Re^wdz 

KiT  ^ 

s^wdz  «  Re^wdz  +  Re^wdz. 


Re^wdz  «  Re^wdz  +  Re\ 
C  S  e 


The  contour  integral  S  follows  a  closed  path  over  that  part  of  the 
slit  on  the  x>axis  which  corresponds  to  the  cavity  hut  excludes  the 
point  z  B  Upon  using  the  fact  that  Re  w  e  u  and  u  b  U^^/2  +  e 
on  the  cavity,  it  is  found  that 


s^wdz  B^udx  B  - y' e  dx  +  « 


RetUvdz  B^udx 
S  S  1 


dx  B  -2c(i  -  1). 


I 
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Furthermore,  from  Appendix  B, 


one  sees  that 


Re  J  wdz  ■  0. 
e 


Using  these  results,  equation  (3«19)  can  be  put  Into  the  form 

Oj_ . 2(2  +  z) (f) ^  fa(/r*^)].^| 


(3.20) 


The  moment  coefficient  about  the  nose  of  the  wedge  Is  defined  as 

C  «  L«y.Bt  to  L 

•  (CHCRD)^  * 


positive  In  the  counter-clockwise  direction.  The  contribution  due 
directly  to  pressure  forces  perpendicular  to  the  x-axls  Is 

<W 

w 


irtiUe  the  y-axls  contribution  Is 


lylay-  c?^ 


C  xdx. 
P 


W 
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&nce,  the  total  moment  C^.  Is 


Cmo  *=  (1  +  xdx, 

W 

hut  since  is  of  second  order  it  is  properly  neglected  in  this 

analysis.  Again,  using  the  results  in  Table  4,  one  has 

CjjQ  =  -(2  +  2)^^  xdx  +  (2  +  (+  ^  +  S/2jxdx.  (3.21) 

W  “  W  '  “  ' 

The  real  part  of  vzdz  is  equal  to  uxdx  on  the  wedge  and 

--fr- 

W  o  1  ' 

Then, 


ReOwzdz  •=  -2llb, 


where  from  eqxiation  (3.9) 


2  I  (T^-l)  (1-T) 


4Tin(>/r +>/^)  \2Mriy  (1+Tj 


As  before,  one  may  write 


Re^wzdz  =  -2llbg  -  Re^wzdz. 


From  Appendix  B, 


►  wzdz  *=  0, 


and  by  using  the  cavity  boundary  conditions, 

/  I  1 


Re^wzdz  K^uxdx  *  "  J  “  l). 

S  S  1  I 

Thus,  Rewwzdz  «  -  -  l). 
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After  introducing  the  above  into  equation  (3>2l),  the  moment  coefficient 
can  be  sinpllfled  to 


0^0  -  [l  . 


(3.22) 


Finally,  one  must  calculate  the  value  of  the  parameter  '^/U^  and 
Introduce  it  into  the  basic  results.  From  Table 


€ 

r 

00 


In  Section  2.3i  it  vas  noted  that  the  quantity  in  the  curly  brackets 
Is  the  cavity  area  A^.  Substituting  the  value  of  A^  from  equation 
(3«llb)  Into  the  above  equation  yields 


«  6  -  1) 

r  '  u/21^.  1)  • 


(3.23) 


Equation  (3 >23)  la  now  introduced  into  the  previously  obtained  results. 


3.1>3.  Discussion 

The  results  of  the  linearized  analysis,  vhlch  eure  summarized  in 
Table  6,  depend  on  the  Independent  variables  £(or  i),  a,  and  relative 
vortlclty  Equation  (3*B)  gives  L  explicitly  as  a  function  of 

Jt  for  fixed  a.  numerical  results,  together  vlth  computation  programs, 
are  tabulated  in  Appendix  F.  Certain  portions  of  these  results  have 
been  plotted  to  Illustrate  the  theory. 

As  predicted  earlier,  the  cavity  length  •  cavitation  numiber 
relation,  the  cavity  area,  and  the  drag  coefficient  are  Independent  of 
the  relative  vortlclty  s/u^.  (ki  the  other  haxid,  the  lift  and  moment 
coefficients  depend  linearly  on  ^/U^j  both  coefficients  are,  of  cotirse, 
zero  in  the  Irrotatlonal,  parallel  flow  past  a  vedge.  For  reference. 
Figures  7  end  8  show  the  drag  coefficient  and  cavity  area  A^  as 
functions  of  Z,  even  though  they  are  the  same  as  those  found  In  an 
irrotatlonal  flow.  The  S  •  Z  relationship  is  plotted  on  Figures  I6 
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and  17;  the  cvirve  for  aajanetrlc  shear  flov  Is  the  curve  labelled 
~  SB  0,  l.e.j  the  Irrotatlonal  flov.  For  this  esse,  Wu  [l6]  has  shown 
tSat  in  the  linearized  theozy,  i  is  limited  to 


Those  results  which  are  directly  affected  by  vortlcity  are  plotted 
in  Figures  9  through  1^.  The  first  figure  shows  e/e  as  a  function 
of  £.  The  next  Is  a  plot  of  a  typical  cavity  shape  which  shows  a 
definite  dependence  on  the  relative  vortlcity.  The  two  vortlcity 
terms  in  equations  (3*13)  acco\nit  for  the  airfoil  shape  of  the  cavity. 
The  second,  or  logarithmic,  term  becomes  large  only  near  the  end  of 
the  cavity  and  tends  to  pull  the  cavity  end  dovnward. 

The  key  results  of  the  theory  are  the  pressure,  lift,  and  mixaent 
coefficients.  Since  the  latter  two  are  linear  In  s/u  ,  Figures  11  and 
12  are  plotted  with  Cj^/(€/U^)  and  functions  of  E 

for  varlovLs  a's.  The  Increase  of  both  coefficients  with  (a)  decreasing 
cavitation  nxmiber,  (b)  Increasing  wedge  angle,  and  (c)  Increasing 
relative  vortlcity  Is  clearly  seen.  As  Z  approaches  zero,  the 
cavity  becomes  Infinitely  long  and  the  magnitudes  of  the  lift  and 
moment  coefficients  approach  Infinity.  Tslen  [12]  found  that  this 
behavior  also  occiirs  In  shear  flow  about  an  Infinitely  long,  solid 
body.  The  pressure  coefficient  Is  presented  In  Figures  13  end  l4. 
The  positive  lift  found  above  Is  represented  here  by  the  area  between 
corresponding  curves.  On  the  xjpper  wedge  surface  near  tte  nose,  the 
pressure  coefficient  of  the  linearized  flov  exhibits  a  large  negative 
value.  This  phenomenon  Is  associated  with  the  hlc^  velocities  required 
for  the  fluid  to  t\xrn  about  the  sharp  nose  point  In  the  equivalent 
nonlinearlzed  flow;  here,  the  stagnation  point  on  the  wedge  Is  below 
the  x-axis  and  behind  the  nose.  In  a  real  flow  a  small  cavity  may 
occur  on  the  \g>por  surface  at  the  nose.  Such  a  cavity  has  been  observed 
in  experiments  with  wedges  at  a  small  angle  of  Incidence  [7,  Chap.  12, 
Ft.  2,  Fig.  12  11.13,  P>33]  where  the  flow  patterns  are  essentially 
equivalent  to  those  of  the  present  rotational  flov.  Figure  1^,  the 


I 
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final  figure  in  the  series,  shows  the  distance  x  of  the  center  of 
lift  from  the  nose  of  the  wedge. 

3*2.  SyianBtric  Flow  Past  a  Wedge 

The  undistvirhed  velocity  profile  of  the  siqpereavitatlng  flow  is 
shown  in  Figure  1  (dotted  line  profile)  and  again  in  Figure  U.  The 
notation,  which  was  Introduced  in  Section  2,  is  shown  on  Figures  4  and  5- 
Note  that  Poisson's  equation  *  s  holds  for  y  >  0,  hut  ■  -e 
holds  for  y  <  0,  contrary  to  the  asynnetrlc  case.  This  change  is,  of 
course,  due  to  the  discontinuity  in  the  vortlclty  at  y  »  0. 

3.2.1.  Solution  of  the  hovindary  value  problem 

It  is  necessary  to  reformulate  the  wedge  boundary  value  problem 
outlined  in  Table  1  because  of  the  sjmmetric  shear  velocity  profile 
at  Infinity.  However,  the  same  conformal  mappings  as  were  used  in 
Section  3.1  may  be  used  here. 

In  Section  2.1,  the  rotational  flow  was  reduced  to  a  harmonic 

2 

or  Irrotationed  flow  by  introducing  the  stream  function  i|rp  ■>  Cy  /2. 

In  the  ease  of  the  symmetric  flow,  it  is  necessary  to  use  the  function 

*W  ■  V*/®’  y  ^  ° 

and  the  function 

tpL  “  y  <  0. 

Then,  one  has  as  desired 

a.  for  y  >  0;  »  «,  Up  »  -«y 

b.  for  y  <  0:  ■  «,  Up  ■  +ey. 

The  remainder  of  the  development  in  Section  2  is  unchanged. 

By  comparing  Figures  2  and  3a  to  Figure  4,  one  can  see  that  the 
boundary  conditions  in  the  synsetrlc  flew  are  (a)  the  saaie  on  the 
upper  and  lower  surfaces  and  (b)  the  same  as  those  for  the  upper  surface 
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of  the  wedge  In  the  asynmetrlc  flow.  The  new  boundary  conditions  on 
the  silt  and  mapped  planes  are  shown  In  Figure  The  mapping  (3*1) 
gives  z  as  a  function  of  ^  as  before.  The  boundary  value  problem 
for  symmetric  flow  in  terms  of  the  con5>lex  perturbation  velocity 
w  =  u-lv  is  the  same  as  that  given  for  the  wedge  flow  In  Table  1, 
except  that  condition  (b)  becomes  Reed,  (w)  ■  U^£/2  +  6,  ^  ^ 
y  =  O".  Note  that  In  the  symmetric  flow, 


Z  «  CT 


It  Is  seen  from  equation  (3 >3)  that  the  function 


w(C) 


2oaj  .  ^  , 

in  £  +  iA(C  -  l/C)  +  +  * 


(3.2it) 


satisfies  the  differential  equation  and  conditions  (a),  (b),  (c),  (d), 
(g),  and  (h)  of  the  revised  bo\indary  ved.ue  problem.  As  before, 
the  remaining  conditions  will  serve  to  determine  A  and  a  relation 
between  i  and  a  for  fixed  €  and  a. 


3.2.2.  Results 

The  function  w(0  can  be  expanded  In  descending  powers  of  z  as 
z  The  series  has  the  form  shown  in  equation  (3>^);  the  boundary 

conditions  (e)  and  (f)  require  that  a^  *  b^  «  a^^  »  0.  By  using 
equations  (3*^)  and  (3*6)  of  Section  3>1>2,  one  has 


w(z) 


OU 
_ c 

“T 


in 


^/^l 

«/C-i 


-A(i-l) 


+ 


1 

“5 

z 


“"c 

.TT 


(M)/r-  ^  (^.1)  (3*h)n«']+  0^^^ . 


(3.25) 


It  follows  that 


A 


I 


-  IK)  - 


and 


OAT 

e 

“  TT 


.  2fJl 


\  V 


+  € 


0. 


IsBuediately,  the  relation  hetveen  a,  tr,  e,  and  i  Is  found  to  he 


^■1  [** 


ViVi  ^  2if'1  Ji 


(3.26) 


The  effect  of  positive  vortlcity  is  to  shorten  the  cavity  for  fixed 
The  complex  velocity  v(a)  can  now  he  vrltten  as 


,  /crtJ  i3/2v 

’'W  ■- 7(^11-) 


(3.27) 


The  remaining  calculations  are  based  on  eq^iatlons  (3.24)  and  (3.27). 
Note  that  in  this  case  the  cavity  remains  aynsetrlc  in  shape  and, 
since  the  flow  Is  symnetrlc,  no  lift  or  moments  can  he  expected.  The 
quantities  of  Interest,  then,  are  the  cavity  shape,  cavity  area,  and 
the  pressure  and  drag  coefficients.  Because  the  flow  Is  symmetric. 


ye(*) 


only  the  upper  wedge  surface  and  \ipper  cavity  shape 
he  calculated. 

The  cavity  area  and  shape  are  found  first.  From  Section  2.3, 
the  hody-cavlty  area  Is 


need 


S  -  +  a  -  -  ^  w(z)EdZ.  (2.1^) 

®  Ti  ■  n 

The  function  is  the  cavity  area  and  a  Is  the  area  of  the  wedge. 
Using  the  contour  paths  In  Figure  6  and  the  analytlelty  of  w(z)  off 
the  silt  leads  to 

^wid*  ■■^w(z)sdc. 

B4C  T 

Fnmi  equation  (3*27)  and  the  theory  of  residues,  one  has 

Ls^wadz  «  - 

OTV 
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Thus, 

and 


S  B 

=  a(i3/2  _ 


(3.28a) 

(3.28b) 


Since  a  positive  vorticity  shortens  the  cavity  length  A  for  a  fixed 
a,  the  area  A  is  reduced  also.  The  cavity  shape  is  easily  found 
by  use  of  equation  (2.lUa)  and  several  results  irom  Section  3.1. 

From  equation  (3.2U),  the  imaginary  part  of  w(0  is 


Im 


Itoffl 


tan“l  C  +  A(C  -  1/C), 


when  C  is  real  (on  the  upper  cavity  stirface).  From  equation  (2.lUa), 

yc  BO  - 

1  ° 

Introducing  the  above  into  this  result,  together  with  equation  (3.12), 
gives 

[?  -  7  -  '•(^1)  t] 

(3.29a) 


As  before,  x  is  given  as  a  function  of  the  parameter  t  by 

.2 


B  i  1  - 


and 


(tST)(tSR) 

T  B  2i  - 1  +  2yi(je-iy 

R  B  2i  -  1  - 


9 


(3.13c) 


Although  the  vorticity  parameter  does  not  enter  explicitly  in 

the  above  equations,  its  effect  is  felt  throu^  the  change  of  i  for 
a  given  a  due  to  vorticity  [see  equation  (3.26)). 

Next,  the  pressure  and  drag  coefficients  are  determined.  The 
pressure  coefficient  is  given  by 

P  ”  P 

C  _  _  ° 


I 
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as  before.  Following  the  same  procedures  as  In  Section  2.3;  one  Is 
lead  to  the  result  that 


Cp  =  (2  +  a)  1^  X  +  <r/2  -  u/U^J  ,  0  <  x  <  1.  (3-30) 

Since  In  this  ease  Is  sysanetrlc  with  respect  to  the  x-axls^ 
eqviatlon  (3«30)  1b  valid  on  both  upper  and  lover  wedge  surfaces. 

From  the  similar  calcxilatlons  in  Section  3*1>2,  It  follows  Immediately 
that  on  the  upper  wedge  surface, 

and 


X  a  i  cos^e/(i-sln^e)  (3*31) 

for  0  <  0  <  Jl/2.  Thus,  equation  (3*30)  can  be  written  In  the  form 
Cj,  -  (a-w)  (x-i)  *  §  (aw)  *  Id  , 

0  <  0  <  n/2.  (3-32) 

The  drag  coefficient  for  a  wedge  Is  given  In  Table  ^  as 

S  -  - 

V 

Utilizing  the  relationships  developed  In  Section  3*1«2,  one  has,  after 
allowing  for  the  present  syioBetry  of  the  boundary  conditions, 

S  ■  ■  35  **  *  (2«)^(-  ^  +  f)  «»]  • 

It  follows  Immediately  from  the  analytic  character  of  v  and  equation 
(3.27)  that 

Da^w^dz  ■  Ai^v^dz, 

W  C 
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since  by  the  theory  of  residues 


"dz  ■  0. 


The  Integral  over  the  cavity  Is  given  by 


Ii^w^dz  =  +  U„a/^U^dy  +y'(6|y^|  +  U^(r/2)yU^dy  I  +  Ii^w^dz, 

C  is  S  J  e 

where  Ln  b  -2uv  on  the  silt .  On  the  cavity  « I y  I  »  and  on 

C 

the  wedge  "ix  =  ^ly^l  =  ®lyl/a.  Thus, 

^lygldy  «  "i^xdy  »  lyidy  *  0, 

W  W  W 

and 

^^ly^jl^y  B€^dy  B  0. 

Furthermore,  from  the  cavity  closure  requirement,  ^  dy  b  o,  with 

W4C 

the  contour  Integral  W+C  being  taken  over  the  wedge-cavity  conibinatlon. 
Hence,  one  obtains 

S  •  ^ 

00  e 

It  Is  again  necessary  to  expand  w(z)  about  z  »  i.  Comparing 
equations  (3.3a)  and  (3.2Ua),  one  sees  that  equation  (B.l)  In  Appendix  B 


becomes 


, .  V  -  2kJ(n) 

w(z) 

vz-/ 


fVom  this  result.  It  follows  that 


‘dz  B 


I 
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Using  the  value  of  found  In  Appendix  B  plus  the  valxie  of  A  found 
previously,  one  has 


C 


D 


(3.33) 


This  Is  formally  the  sane  equation  fotind  In  the  Irrotatlonal  case) 

however  the  cavity  length  i  in  equation  (3*33)  Is  altered  by  a  first 

order  vortlclty  effect.  This  alteration  of  £  causes  a  displacement 

of  the  curves  of  Cjj  versus  a. 

In  order  to  define  a  value  for  €  to  complete  the  solution  of 

the  symmetric  problem,  one  must  resort  to  physical  and  mathematical 

Intuition.  Since  the  circulation  P  Is  Identically  zero  about  the 

cavity  In  the  symmetric  flow,  one  cannot  \ise  the  method  of  matching 

circulations  which  was  so  successful  In  Section  3*1«2.  However, 

since  c  Is  to  be  a  representative  valiae  of  cjy  I  over  the  whole 

c 

cavity.  It  Is  reasonable  to  choose 


l.e.,  the  mathematical  average  value. 


3.2.3.  Discussion 

The  results  of  the  symmetric  flow  analysis  are  summarised  In 
Table  7*  Appendix  F  lists  the  numerical  computations  which  have  been 
carried  out  by  using  the  equations  In  Table  7*  Three  important  results 
are  presented  In  graphical  foxm  in  Figures  16  through  I9. 

First,  the  cavitation  number  a  Is  plotted  as  a  function  of 
cavity  length  i  In  Figures  I6  and  17.  It  Is  evident  that  a  positive 
vortlclty  €  causes  a  induction  In  the  cavity  length  for  fixed  values 
of  O'.  Although  the  cavity  length  Is  Infinite  ^en  a  0  In  an 
lrrotatl<mal  flow  (c  >  0),  there  Is  a  maximum  I  corresponding  to 
each  value  of  the  relative  vortlclty  t/u^  >0  in  a  rotational  flow. 
When  c/u^  <  0,  the  solution  to  the  problem  Is  no  longer  unlqvie. 

Then,  as  seen  In  Figures  I6  and  I7  or  In  the  tabulated  data,  there 
exists  a  minimum  a  for  each  value  of  0.  Fbr  each  a  greater 


than  the  mtniimim  value^  there  are  tvo  posslhle  cavity  lengths  -  the 
conjugate  lengths.  When  c  Is  less  than  the  mlnlnum  value,  no 
solutions  exist  to  the  linearized  problem.  When  tvo  solutions  exist, 
both  satisfy  all  imposed  boundary  conditions  and  produce  physically 
reasonable  drag  coefficients  and  cavity  shapes  (see  Figures  l8  and  19). 
In  spite  of  the  large  vortlclty  effects,  the  lover  limit  on  the  cavity 
length  still  seems  to  be  that  value  determined  for  Irrotatlonal  flov 
(see  Section  3*1»3)  and  denoted  by  1  ^ 

Second,  the  drag  coefficient  Cjj  Is  given  as  a  function  of  c 
and  oe  for  tvo  values  of  In  Figure  l8.  A  comparison  of  this 

figure  vlth  Figure  7  shovs  that  vhen  the  vortlclty  Is  positive,  the 
values  of  for  fixed  a  and  o  are  slightly  Increased  over  the 
corresponding  Irrotatlonal  valties.  Hovever,  vhen  €  <  0,  the  drag 
coefficient  Is  reduced  and  the  Cjj  curves  In  Figure  l8  lie  below 
the  comparable  cvurves  In  Figure  ?•  One  should  note  the  reappearance 
of  a  minimum  cavitation  number  for  €  <  0  and  the  tvo  i>osslble 
valijes  for  each  c  above  the  minimum. 

Finally,  Figure  19  shovs  the  shapes  of  the  cavities  trailing 
behind  a  vedge  vhen  the  cavitation  number  Is  fixed  and  c/U^  Is 
varied.  The  tvo  longest  cavities  shown  are  the  conjvigate  length 

cavities  for  «/U  =  -  0.080. 

•0 


rlc  Flow  Past  a 


The  final  problem  considered  Is  a  parallel,  uniform  shear  flow 


past  a  supercavltatlng  flat-plate  hydrofoil.  The  unit  length  hydrofoil 
Is  placed  at  an  angle  a  with  the  x-axls  (as  shown  In  Figure  l). 

The  developnent,  based  on  the  methods  outlined  In  Section  2,  closely 
follows  Parkin's  solution  of  linearized  cavity  flow  past  a  hydrofoil  In 
a  liquid  with  gravity  ].  The  similarity  between  the  present  case 
and  the  gravity  flow  problem  will  be  evident.  In  fact,  the  basic 
boundary  value  problems  for  the  tvo  flows  differ  only  by  a  change  In 
sign  of  the  pertxizbatlon  paraiseter,  l.e.,  </U^  ■  -  where  u/U^ 

represents  Parkin' s  gravity  pazasieter. 


I 
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3*3*1*  Solution  of  the  boundary  value  problem 

As  before,  the  confomal  napping  technique  la  used  to  solve  the 
bouxidaxy  value  problem  given  In  Table  1.  The  transformations  are 
listed  In  Table  2;  In  Figure  3b  the  silt  s-plane  and  transformed  planes, 
together  vlth  cozrespondlxig  boundary  conditions,  are  shown. 

The  solution  v(C)  Is  constructed  from  the  singularities  In 
Table  3.  By  conqparlng  the  form  of  the  boimdary  conditions  axid  the 
available  singularities,  one  can  see  that  v(C)  should  have  the  form 

v(C)  -  tA  H)  +1B  +  iCin  C  +  1^  +  E.  (3.35) 

The  boundary  conditions  In  Table  1  are  applied  to  equation  (3*35)  to 
determine  the  real  constants.  suBaaary,  this  leads  to  the  following: 

a.  Re(v)  >  U  £/2  +  6  for  C  and  >  1,  so  that 

E  «  U^Z/2  +  6. 

b.  Re  (w)  «  U^E/2  -  6  for  C  and  5  " 

-  HC+E  =  U^E/2  -  e  and  C  «  2e/ll. 

c.  3hi(v)  «  -KAJ  for  C  on  bhe  unit  semi-circle,  so  that 

c 

-  ^  +  B  «  CAJ^  axvl  B  »  CAJ^  + 

d.  v(b)  0  as  z  ■*  -w.  In  the  (^plane,  as  z 
Q  Ik  with  k  «  lair,  as  before. 

Equatlw  (3*35)  becomes 

w(C)  -  iA(C-l/C)  +  iou^  +  ^  Ax  C 

+  ID  (i  +  +  U^2/2  +  c.  (3.35a) 
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In  this  form,  v  satisfies  the  differential  equation  and  all  boundary 
conditions  except  (e)  and  (f)  in  Table  1.  These  remaining  boundary 
conditions  allow  determination  of  A  and  D  and  development  of 
a  cavitation- number  -  cavity- length  relationship. 

Following  Parkin  [U],  one  couqtletes  the  solution  of  the  problem 
in  the  Q- plane  (see  Figure  3b).  From  condition  (f)  and  item  d  above, 
w(z)  ->0  as  z  -*■  -»,  while  in  the  Q-plane,  Q  -►  ik  as  z  ■*  -«o. 

Equation  (3* 35a)  may  then  be  written  in  terms  of  Q,  and  Q  must  be 
allowed  to  approach  Ik,  where  w(Q)  approaches  0.  La  accordance 
with  the  transformations  of  Table  2, 

4Q  +  2  »  C  +  1/^; 

hence, 

4  =  1  +  2Q  +  2V4(ftfl)'.  (3.36) 

The  negative  root  is  chosen  for  Q  real  and  <  -  1;  otherwise,  the 
positive  root  is  appropriate.  When  Q  «  ik, 

C  »  2ik  +  1  +  2yikIF. 

Letting  r  »  euad  s  «  one  obtains  from 

equation  (3 •35a)  and  the  above 

w(lk)  «iAri+yk'8  +  i(2k  +  y? r)  -  ■  ^  1  +  icflj 

L  l+ZiTs  +  l(2k+inrr)J  ® 

+  iD(i  +  [i^s+11^rr^)+l(2€/n)ln((l+yk’s)®  +  (2k  +yk'r)^] 

-  (2€/n)  tan‘^  +  U  E/2  +  e  -  0.  (3.37) 

1  +yi?  s 

lotting  A  B  A  U  and  D  ■  DU  as  Parkin  does,  one  obtains  two 
equations  for  A^  and  D^  fr<x&  the  real  and  imaginary  parts  of  equation 
(3>37)>  These  are 
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The  equations  are  solved  simultaneously  to  yield 


‘  ik' t  {■  t  ' 

i}] 

(3.38a) 

and 

-f}-! 

fe  1  . 

l^d+yiTs) 

+  (2k  +  Vk*r)^J  +  • 

(3.38b) 

In  the  Q-plane,  using  equation  (3.36),  the  conplex  perturbation  velocity 
Is 


v(Q)/U^  -  i4A^Vft(Qfl)'  +  iD^1ll[Qfl)/Q/2  +  la(l+Z/2)+E/2 

+  i/u^  +  1  J-  I  Aa  [l+2(Jf2lfe(Qfl)’l.  (3-39) 

•0 

The  remaining  condition  (e)  Is  the  closure  condition.  As  noted 
previously,  if  the  cavity  Is  to  close,  the  net  strength  of  sources 
vlthin  the  cavity  sust  equal  sero.  It  is  equivalent  to  require  tiiat 
v(z)  have  no  real  residue  vithln  an  infinitely  large  circle  T 
(see  Figure  6)  surrounding  the  cavity,  i.e., 

&i^v(E)ds  >  0. 

T 
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Because  v(z)  Is  harmonic  off  the  slit  in  the  z- plane,  it  follows  that 


Ln 
B+C 


^v(z)dz 


0, 


(3.**0) 


with  B+C  the  boundary  of  the  foil-cavity  system.  Thus,  in  the  Q-plane, 
the  closvire  integral  is 


Ic  =^w(Q)  ||  dQ.  (3.41) 

B+C 

The  foil-cavity  system  extends  from  -«  to  +<»  in  the  Q-plane.  This 
Integral  is  evaluated  in  the  classical  manner  by  using  a  semi-circle 
of  radius  R  in  the  upper  half  plane  to  form  a  closed  contour 
One  has,  from  the  given  transforms. 


dz 

d5 


2ik^(k^  +  Q^)“^Q« 


(3.42) 


Thus, 


and 


R 


dQ  -►0 


as  R 


I 

c 


2111  (Residues  within  Cp).  (3.4la) 


The  minus  signs  account  for  the  reversal  of  the  line  integral  orienta¬ 
tion  in  the  Q-plane.  The  only  residue  within  Cp  occurs  at  the 
second  order  pole  Q  »  Ik.  The  residue  at  the  pole  Is  given  by  [21] 


b 


1 


ik  ^ 


(3.43) 


since  w(ik)  >  0.  The  introduction  of  equation  (3.39)  i&bo  equation 
(3.43)  and  the  use  of  the  subseqwnt  result  give 


irfi 


s(l6k  A  +D^) 
0  0 


-  8kr^A^4  y  "  ^  ( V  i  * 


(3.44) 
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From  equations  (3«Ho)  and  (3. 1*1),  It  follows  that 


-  r 


8k8(A^  +  J  ^  )  +  r(l6k2A  +  D  )  «  0 
«0  / 

for  closiare.  By  Introducing  the  values  of  the  constants  A^  and 
one  derives 


in [(1+1^  s)^  +  (2k+1^  r)^]  +  a  (l+^ 
fi  /,  2  ^  .-1  2k+lfl?r\  .  zl  V  Sksi  -i/I:/  fi  T  2  ^...-l  2k+VPr“l 

^  ^  Ul*VPef  +  (2kHVrr)^]  +  a(l+||Jj  «  0. 


(3.45) 


3. 3*2.  Besults 

the  prevloxis  Section  3*3*1>  the  basic  solution  v(Q)  and  a 
relation  (equation  (3*^^))  between  the  problem  parameters 
i,  £,  and  a  were  found.  Based  on  this  information,  the  cavity 
characteristics  —  length  (closure  condition),  area,  and  shape  — 
may  be  found.  In  turn,  one  can  determine  the  pressxire  force  coefficients 
azkd,  finally,  the  vorticlty  parameter 

I^n  simplification,  equation  (3>^^)  yields  the  closure  condition 


a 


k+yrj; 


(3.46) 


As  C^/U^)  '*  0,  this  equation  reduces  to  Tulin's  condition  for  clostire 
in  the  Irrotatlonal  case  (6),  l.e., 


a 


k£ 


I 
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with  S  *  a  when  "i/u^  s  o.  When  "i/U^  is  not  zero,  a  positive  e/U^ 
produces  a  lengthening  of  the  cavity  over  the  irrotational  case  for 
fixed  a  and  Z. 

The  cavity  shape  is  found  by  integrating  eqviation  (2.l4a)  from 
the  appropriate  point  on  the  hydrofoil,  i.e.,  the  leeuling  or  trailing 
edge.  By  enqploying  the  previously  determined  relation  for  dz/dQ  and 
the  fact  that  z  and  Q  are  real  on  the  cavity,  one  can  Integrate 
equation  (2.l4a)  directly  in  the  Q- plane.  For  real  Q, 

Im  ^  +  D +  a(^)+  in  |l+  2Q  +  2VQ^ftfl)  |  , 

00  '  '  00 

on  the  upper  and  lower  cavity  surfaces  respectively.  It  follows  that 

a.  on  the  upper  cavity  surface,  for  q  >  0  and  real, 

^f\_  +  r  ^  (2+Z) 

+  (iin|l+2Q  +  2VQ(Q5)jj  .(3.47a) 

b.  on  the  lower  surface,  for  q  >  1  and  real, 

1  -  ^lfe(Q-lJ  -  ^  (2+2) 

.^Qin|l-2Q.2l^|]^-5^. 

(3.U7b) 

By  use  of  the  tabulated  integrals  from  Appendix  C,  these  equations 
may  be  written  as 


^c 

i 


(3.48a) 


on  the  upper  surface  and 


a 

r  "’7  +  ?+2 


(“•Vi.  -  ^  - 


%•  I,  . 

^  ^  l*i(k^  +  q2) 


(3.1^8b) 


on  the  lover  cavity  surface.  The  value  of  x  for  a  corresponding  q 
Is  found  by  Integrating  equation  (3 >^2)  and  Is  given  by 

X  .  (3.J+9) 

7  “  — j>  '  6 

^  (k2  +  q2) 

Since  the  cavity  and  body-cavity  areas  are  the  same,  the  cavity 

area  A  Is  [from  equation  (2.1^a)] 
c 


‘  1+z 


k  A 


The  contour  Integral  Is  again  evalviated  In  the  Q-plane,  where 


^2^ 


adz  a>  -2k 


As  before,  the  theory  of  iresldues  Is  applied,  nie  only  pole  In  the 
region  Is  the  third-order  pole  at  Q  ■  Ik;  the  residue  there  Is  given 
by  [21] 

^  k^/^  d^  rw(Q)Q3n 


^  d?  LcOflk)?]  ^  ^  ‘ 

This  equation  Is  reduced  In  Appendix  D.  The  Imaginary  part  of  the 
result  Is 

^c  “  +  (2k^  +1)8]  -  ^  ‘  TO  * 

(3.50) 

with  A^  and  given  by  equations  (3.30). 

The  calciilatlon  of  pressure  force  coefficients  depends  on  the 
evaluation  of  the  pressure  coefficient,  which  represents  the  difference 
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'between  the  pressure  p  on  the  lover  surface  of  the  hydrofoil  and 

the  cavity  pressvure  p  on  the  upper  surface.  These  relationships 

c 

were  defined  In  Section  2.3  and  tabulated  In  Ta'ble  4.  From  that  ta'ble, 

0p  =  -(2+Z){f  (3.51) 

In  the  ^plane  the  lower  side  of  the  foil  Is  represented  "by  Q  *= 


0  <  I  <  1*  Also,  u/U^  «  Real  (w/u^).  Using  the  above  relations 
and  equation  (3>39)>  has 


^  =  -  4AqV|(1-0'  +  ^ 


26 

s  ■  W 


since  z  e  iQ^/(k^+Q^),  then  x  «»  on  the  hydrofoil,  and 

C  can  be  written  as 
P 


-  (2.1) 


] 


iC  2 

2  2  "  1 
k  +1 


0  <  5  < 


(3.52) 


with  tan”^(  )  <  n. 

The  normal  force  coefficient  Cj^  Is  given  directly  by  Integration 

of  C  over  the  foil.  Thus,  from  Table  4, 

P 


'N 


o 


and,  from  equation  (3<^l), 

(3.53) 

'  "  O 

The  Integral  In  this  equation  can  be  written  In  terms  of  a  contour 
Integral  aroxmd  the  hydrofoil-cavity  system,  l.e. 

I 


A  coiq>arlson  of  the  velocity  functions  used  in  the  hydrofoil  and  wedge 
flows,  together  with  the  results  of  Appendix  B,  show  that  there  is  no 
contribution  to  the  normal  force  from  the  Integral  of  w  around  the 
end  of  the  cavity.  This  is  also  true  for  the  integral  wz  about  the 
end  of  the  cavity.  Since  u/U  «  Z/2  +  "i/U  on  the  upper  and  lower 
cavity  surfaces  respectively,  one  obtains 


The  introduction  of  this  result  and  equation  into  eqviatlon  (3>^3) 

produces 

Cjj  -  -(2+2)  +  Re(lg/U^)J  . 


The  final  form  for  is  obtained  by  substituting  the  value  of 

from  equation  (3*^^)  and  the  values  of  the  constants  A  and  into 

o  0 

this  result.  After  some  manipulation,  one  has 


2+2) (fCk)  ^Qk(2+2)  +  Z  +  in[(l-»^s)^+(2k+l/Fr)^] 


4.  1 


Where  from  Table  1»,  ■  Cjl>  ■  oeC^,  and  D/L  ■  a. 

The  moment  coefficient  about  the  leading  edge  of  the  hydrofoil 

is,  from  Table  4, 


S«0 


o 


When  C  is  Introduced  into  this  result,  one  derives 
P 


‘'MO 


-(2+2)  -  I  +  xdxj  . 


(3.55) 
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As  before,  the  Integral  on  the  foil  Is  transformed  into  a  contour  inte¬ 
gral  with  the  result  that 


B+C  " 


The  integral  in  this  eq\iation  has  been  evaluated  in  Appendix  D  in 
connection  with  determination  of  the  cavity  area  A  .  Taking  the 
real  part  of  that  restilt  and  introducing  it  into  the  above  equation 
prodiK:e8 


^MO  (2+2)^-  ^  |iA^([6i-l]r-[l2i+l]ks) 

■^2^  +  £1  +  ITT  (t**k^+5lr+ks)J  J  . 

After  minor  re-arrangement,  the  moment  coefficient  has  the  form 


“'MO 


=  (2+2)  (t41fl]r+ks)-i^  +  ^ 

*  •  ([6/-l]r.[12*n)k.)+  ^  ^r-  [»  +  i]  • 

(3.56) 

Finally  one  determines  the  vorticlty  parameter  l/U  .  From  Table 


(3.57) 


It  has  been  noted  that  this  result  can  be  written  directly  as 

^  e(A^-a/2)  6(A^-a/2) 

00 

mhere  k  e  Vi-1.  Recalling  equation  (3*50)  which  gives  A  ,  one  can 

C 

see  that  e/u  may  be  found  directly.  This  calculation  involves 


I 
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solving  a  quadratic  relation  obtained,  by  cosiblnlng  equations  (3>^) 
and  (3.^7)  and  using  equation  (3>U6),  the  closure  condition,  to 
eliminate  £.  The  result  obtained  Is 


U 


(i) 


i— 4 - 

*  -  -1 

r-kXi,+(k-a)X.,  Ip  4ca.  M 

yL  r  ^  ^  j 

(3.58) 

in  terras  of  the  notation  introduced  in  Appendix  E.  For  €  <  0,  the 
minus  sign  is  taken  to  preserve  the  form  of  */U^*  This  is  required 
because  <  0  when  €  <  0.  With  ®/U^  known  explicitly  as  a  function 
of  i(orZ),  and  a,  the  solution  Is  complete. 

3.3*3«  Discussion 

The  solution  to  the  hydrofoil  problem  Is  summarized  In  Table  8. 
Appendix  F  lists  a  sample  of  the  numerical  computations  carried  out 
by  using  the  equations  In  that  table.  The  resialts  are  presented 
graphically  In  Figures  20  through  30. 

The  cavity  length  i  is  plotted  as  a  function  of  cavitation 
number  £  for  two  values  of  the  attack  angle  a  In  Figures  20  and 
21.  The  effect  of  vortlcity  In  the  flow  is  clearly  seen;  for  a  given 
L  the  cavity  Is  lengthened  by  positive  vortlcity.  In  Figure  22, 
the  effect  of  vortlcity  on  the  size  of  the  vortlcity  parameter 
Is  shown.  A  typlccd  cavity  shape  Is  plotted  In  Figure  23*  The 
effect  of  positive  vortlcity  Is  em  increase  In  the  cavity  width  and 
the  cavitation  number  for  a  fixed  cavity  length.  In  contrast,  Parkin  [U] 
found  that  cavities  effected  by  a  transverse  gravity  field  lie  Inside 
the  corresponding  gravity-free  cavity. 

The  normal  force  and  moment  coefficients  are  pictured  In  Figures 
24  through  27.  (Recall  that  «  Cj^  and  Cp  ■  ^  coefficients 

are  plotted  versus  £  In  Figures  2k  and  26  for  a  relative  vortlcity 
e/U^  >  0.04.  As  the  cavitation  nuniber  approaches  zero,  both  sets  of 
curves  turn  upward  and  Increase  rapidly.  This  rapid  increase  In  lift 
and  moment  as  the  cavity  becomes  infinite  In  length  Is  consistent  with 
the  results  of  Tslen's  Investigations  of  shear  flows  [l2]  and  the 


I 
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results  of  Section  3*1  (see  Figures  11  and  12  end  discussion  on  Page  37)* 
Figures  2?  and  27  show  and  as  functions  of  c/u^.  Oenerally, 

both  Increase  as  e/u^  Increases;  however,  the  increase  becoiaes 
pxvnounced  only  as  Z  becomes  small.  The  variation  of  the  location 
of  the  center  of  lift  as  a  function  of  £  and  e/u^  is  given  in 
Figure  28. 

Two  sets  of  typical  pressure  coefficients  are  presented  in  Figures 

29  and  30.  The  first  shows  the  effect  of  vortlclty  on  at  constant 

cavity  length.  The  second  shows  the  effect  of  the  cavitation  nvaiber  Z 

on  when  the  angle  of  attack  and  vortlclty  are  held  constant.  In 

this  shear  flow  theory  no  pathological  cases  (l.e.,  cases  in  which 

turns  sharply  downward)  are  found  which  would  coopsure  to  those 

eigiwrienced  by  Parkin's  gravity  theory.  The  reason  for  this  difference 

is  that  the  change  in  sign  of  the  coefficient  D  (Parkin's  A  ),  which 

o  o 

causes  negative  lift  in  the  gravity  case,  doesn't  occur  in  the  present 
problem  because  the  perturbation  parameter  S/U  differs  from  Parkin'  s 
gravity  parsawter  u/U^  &  minus  Biffi,  Thus,  while  a  strong  gravity 
Influence  causes  negative  lift,  a  strong  vortlclty  Increases  the  lift. 

Finally,  from  Figures  20  and  21  and  the  tabulated  data  it  is  seen 
that  when  c/u^  >  0,  the  problem  solution  is  not  uniqtue.  Ttere  is  a 
mlnlsnai  Z  for  each  value  of  ^  0.  When  Z  is  greater  than  the 

minlaium  value,  there  are  two  possible  cavity  lengths-  conjugate  lengths; 
when  Z  is  less  than  the  miniiBua  value,  no  solutions  exist  to  the 
problem.  As  in  the  ease  of  symmetric  flow  past  a  wedge,  when  the  two 
soluticms  do  exist,  they  both  satisfy  all  conditions  of  the  problem 
and  produce  physically  reasonable  joessure  force  coefficients. 
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4.  COMCLUDIHG  REMABKB 


In  Section  3  the  linearized  theory  has  been  applied  to  three 
rotational,  supercavitating  flows.  From  these  applications  it  may 
be  conclvided  that  the  effects  of  rotation  (vorticlty)  are  significant 
whenever  the  magnitude  of  the  relative  vorticlty  e/u^  is  greater 
than  0.02.  Experimentally,  the  vorticlty  effects  are  most  likely 
to  be  detected  in  measurements  of  lift  and  moment.  In  the  cases  of 
symmetric  shear  flow  with  negative  vorticlty  and  uniform  shear  flow 
past  hydrofoils  with  positive  vorticlty,  further  analysis  and  esqperl- 
mentatlon  will  be  required  to  determine  (a)  if  the  non-unique  solutions 
found  in  Sections  3*2  and  3>3  ^  occur  and  (b)  if  the  cavity  ml^t 
tend  to  oscillate  between  the  ccm Jugate  lengths  and  hence  cause  some 
dynamic  effects. 

In  the  uniform  shear  flow  about  wedges  and  flat-plate  hydrofoils, 
the  positive  rotation  has  been  shown  to  cause  an  Increase  in  lift  and 
moment  forces.  In  hydrofoil  flows  an  attendant  Increase  in  the  drag 
can  be  expected.  It  is  also  important  to  recall  the  large  increase 
in  the  size  of  vorticlty  effects  which  occur  as  the  cavity  lengthens. 
The  work  of  Parkin  and  others  on  associated  pzmblems  suggests  that  the 
present  linearized  theory  may  over-estimate  the  vorticlty  effects 
when  the  cavity  is  extremely  long;  however,  the  theory  gives  no 
Indication  of  failure  in  these  regions.  But,  as  the  cavity  length 
approaches  infinity,  the  lift  and  mnsent  coefficients  do  become 
infinite.  On  the  other  hand,  in  the  symnetrlc  flow,  the  cavitation 
nuniber  and  cavity  length  as  well  as  the  effects  of  vorticlty  were 
genezmlly  found  to  be  bounded. 

The  present  application  is  limited  in  several  ways.  The  wedge 
half-angles  and  the  hydrofoil  attack  angles  are  bounded  by  the  \;vper 
limits  associated  with  the  basic  linearised  theory.  Chen  [2]  notes 
that  the  linear  theory  predicts  the  force  coefficients  with  an  error 
of  6  percent  for  a  flat-plate  hydrofoil  at  ^  degrees  incidence  and  an 
error  of  ^  percent  for  a  symmetric  wedge  of  1^  degrees  included  angle. 
Finally,  there  is  the  method  of  determining  the  vorticlty  effect. 

The  method  chosen  is  arbitrary  to  be  sure,  but  its  value  lies  in  the 
fact  that  the  method  (a)  permits  a  eonparitively  simple  solution  of 
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an  othervlse  difficult  problem  cuad  (b)  accounts,  In  general,  for  the 
over-all  vortlclty  effects  rather  than  the  effects  at  one  particular 
point  In  the  flov.  A  specific  objection  which  may  be  raised  Is  that 
the  present  theory  cannot  be  extended  directly  to  the  second  order 
because  of  the  averaging  technique  which  Is  used.  It  Is  anticipated 
that  the  present  method  can  be  refined  In  future  analyses  and  also 
extended  to  more  general  rotational  flows. 


I 
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AFFENDn  A 
Singular  Integral 


The  solution  of  the  singular  Integral  equation  of  the  first  kind 
i 

-  f  (x),  1<  x<  i,  (A-1) 


Is  accomplished  by  reducing  the  Integral  equation  to  one  of  known 
type  by  a  change  of  variables.  Let 


1  and  r  ■ 


S^.i, 


so  that 


Then, 


dt,  dx  «  dr 


I  -  i[(i-l)(t+l)  +2],  X  .  i[(i-l)(r+l)  +2]. 


|i(l)  *►  0  (t)  and  f  (x)  -►  f  (r), 


while  equation  (A-l)  becomes 


(*-2) 


This  equation  occtirs  In  theory  of  airfoil  motion  and  Is  often  called 
tte  "airfoil  equation."  Rote  that  Improper  Integrals  must  be  taken 
vlth  their  Cauchy  principal  value.  The  solution  to  (A-2)  Is  provided 
by  Trlconl  [22];  It  Is 


?  (t)dt .  — L , 


(A-3) 
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where  represents  the  non- trivial  solutions  of  the  homogeneous 

equation 


-1 

The  solution  to  (A-2)  is,  as  expected,  non- unique;  the  constant 
is  determined  by  the  condition  of  smooth  separation  at  the  trailing 
edge  of  the  wedge.  Returning  to  the  original  variables, 

a 


^l(x) 


rvfi-xXx-iy 


I f  (8)^8  ♦  c  ,  (A-k) 


where  C  has  replaced  2S’/ll^(i-l)  to  achieve  the  desired  form. 
For  ^(x)  to  remain  finite  at  x  «  1, 

i 


f  (5)d|. 


Thus,  the  distribution  function  |i  can  be  written  as 


m(x) 


(A-5) 
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APPENDIX  B 
Ebcpaiuiloas  at  z  ■  i 


la  order  to  calcvGLate  the  pressure  force  coefficients  in  a  wedge 
flow.  It  Is  necessary  to  expand  the  complex  perturbation  velocity 
about  z  s  i  and  to  calc\ilate  contour  integrals  over  the  circle  e 
(see  Figure  6).  The  velocity  v(z)  becomes 


w(z) 


^  ^  ^  2Vi(je-l)'  -  iniz^  (B-1) 

Vz  -  i ' 


as  z  ■>  i  because,  In  equation  (3>3&)> 


VTTT 

irf^  -►  ID 

r+1 


a.  C 


b. 


>  Id 


c.  iA(c  +  o(vr7) 

’  Tt-A 

d<  Jkx  0 

Thus,  v(z)  may  be  written  as 


vu)*K.ssjsni. 

irrj  " 

where  from  above  K  a  0(a,  «).  nrom  equation  (B-2) 

IT  ^zTT 


(B-2) 


w2  (.) 


.  jb  VTT*  SJSffiSl  A. 

^  nCT 

(B-3) 

?lntefi7als  whl<^  occur  In  the  linearized  theoz^  are  of  the 
wdz,  wwzdz,  Qw^dht.  Since  only  singular  terms  may  contribute 

to  sucS  contSur  intee^als,  one  needs  to  evaluate  contour  Integrals 
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for  the  following  singular  fimctlons 


b.  fg  «  in  Vz  -  / 

C.  f^  »  z/l/z  -  i 

d.  f ■  z  inVz  “  ^ 


e.  f^  «  (inV  z  -  i)/i^  -  i 

f.  fg  »  (in  Vz  - 


g.  f. 


7  z  - 


^  A  TM 

Near  z  ■  i,  z  -  i  ■  re  and  z  ■  re  +  i,  where  r  -►  0  as  z  -►  i. 
Id 

Then,  dz  »  Ire  dd  on  e  and  an  Integral  is  given  by 

n 

Ji  -  yl^B^  e^®d0. 


.10 


.n 


All  integrals  of  the  fora 


u 

E  e  J  f^(e^®,  0)d0 

-n 


are  finite  and  if 


•'l  •  r‘  E,  t  >  0, 


then  such  ■  0. 

Introducing  the  above  notation  into  the  f^  and  writing  the 
corresponding  gives 


.n 
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il  il 

To.  Jg  =  Ir  in  J'e^®de  -  ee^®  de  «  0 

-n  -n 

since  yl^BQ  r  Aa  ■  0  for  any  finite  t. 

n 

c.  de  +  iJ^  «  0. 


-n  n 


d.  Jj^  =  ir^  in  J e^^®  dS  -  ^.1^  ^Jee^^dB  +  ijg  =  0. 

e.  =  ^1^  ir^/^  in  e^®/^  de  «  0 


since  j.1^Dq  r^^^  in  ■  0. 


f. 


«  ^1^  ir  J|^(in  r^/2)2  +  ie  Ai  r^/^  -  |^J 


ei®de. 


Since  all  other  terms  go  immediately  to  zero,  one  has 


e^®  dfi. 


By  L* Hospital's  Rule 

._  (fa  ^ 


’  r^^  ‘l  “  ^6  “  °* 

n  "  “5 


8-  *^7  *  r^-H 


It  foUovs  that 


^vdz  B  0, 
e 

^  vzdz  =  0, 


«:^w^dz  *  2ni(4A^i(i-l)), 


(B-U) 


I 
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AFFEHDIX  C 

Svnaaary  of  Cavity  Shape  ]&itegralB 


The  results  listed  below  and  the  accoiopanylng  definition  of  terms 
are  susanarlzed  from  reference  [U].  Appropriate  changes  of  notation 
have  been  made  to  make  the  results  cosgpatlble  with  this  development. 
The  new  terms  defined  are  as  follows: 


(a  «  kVk^+l 


03  03 


3c 


■y 

7,  =-5 

03 

One  then  has 
^  ^  (k%Q^) 


V^> 

kPg 

k^ 

^  2ai 

and  ■  —-5  . 

203 

r  2(kW) 


L - 

.  ^  L-1  i .  .  I  ^]\ 

^ dQ  1  /dfoT^TT 

h’J 

— J 

. X  fw" I .  _ - 1 .t».- 


i}' 


I 
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=  /  ^n|l+2ftf2yQ(ftfl)’|QdQ 

in[l+2a+2yq((i+l j  ]  ^  f  ^  p  1 

a(k=-fq®)  a?Y‘ ”[  (?'*~?)(6fT^)  J 

*  ^  L-i  I  *  t„-i  ,  -  i  -t«-^  ^  I) , 

"  \  ’  r^q*8l-tro^L7  ^  *1// 

^  /ajS^TL 

^  (k^  +  Q^) 


2 


,  ■ «  .  1  r  ^  ^  rar,q.v,g?^3^.[.r,q^,l-)/1 

2(k®n®)  *  L  ''  L(qW)(tv6^)2+t.r2+»g)2j  J 

.  ^  (t.n-"  I  .t«.-i  ;i^)l , 

“  \  ■*  V- VvSn^  ^  /J 

■fw^ 


^  1  fqlC?7I)',^i  .  f(t73^q-6i-Vq(q-i)1^+t-rgq*5  ]2)/  "j 

=  .t.n-^;2^j\. 

q-1)  1  1  /J 


-  \  f  t«.-i  i  .t«-i 
2“  \  ’  r^q-V 


and 


afali-aa-afcftt-D'Ida 

^  fall-2Q-2lfa(q.l)l  1  [^8 

2(k^  +  q®)  ^  [(qWxtr^-s^l^+t-rg+Sg]®)  J 


I  -  tqn"*-  k  +  tm"^  * 

"  I  ’  r.q-V 


Vq(?i)’ 


-  tan 


.l-’’2^2V 

■ 


I 
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APPENDIX  D 

Evalvtatlon  of  CoMPlex  Int«gr»l 


In  the  hydrofoil  analysis  of  Section  3*3>  It  Is  necessary  to  find 
the  real  and  Imaginary  parts  of 

zd,  .  .  49.  (D-1) 


cavity 

From  the  theory  of  residues^ 


c  -  -  anu, 

C  1 


vlth 


h^  B  the  residue  at  the  triple  pole  Q  >  Ik.  Nov, 


Ik* 


(D-2) 


On  performing  the  Indicated  differentiation  and  noting  that  v(lk)  ■  0, 
one  obtains 

'’i  ■  ^  (t )  (t )  J 

From  equation  (3*39)^  the  derivatives  of  the  complex  velocity  are 

emd 


^  ■  nfe  [■“»  *  ^  ■  k  ‘“^"0  • 

These  derivatives  are  evaluated  at  Q  >  Ik  and  Introduced  Into  equation 

(D-3)*  resulting  value  of  b^  Is  then  s\d>8tltuted  Into  equation 

(1^2).  After  some  simplification,  Including  the  use  of  equation  (3.44) 

and  the  closure  condition,  C  may  be  written  as 

c 


I 
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^  els’  els' 


^A^[(61-l)p-(l2i+l)lM]  +  ^[p-(2k  +  |)8] 

[(Uk^+5)r+k8]  +  1  ^4A^(k[Uk^+5lr+[2k^+l]8)  -  •gS  ^[2k+|]r-8  J 

(kr-8)j^  . 
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APPEHDn  E 

Solution  of  Quadmtle  Equation  for  c/u 


The  vortlclty  paraneter  is  found  as  the  solution  to  a 

quadratic  equation.  The  equation  la  obtained  by  eoeiblnlng  eqmtlona 
(3.^)  and  (3.^7)  and  using  eqxiatlon  (3*^6),  the  closure  condition, 
to  eliminate  £.  Parkin  has  solved  the  analogovis  gravity  paraneter 
problem  [4].  The  nonenclature  tised  here  Is  essentially  Parkin's, 
but  the  results  are  sonsidiat  different. 

The  follovlng  simplifying  and  systematising  terms  are  Introduced: 


H(k) 

Ci(k) 

C3(k) 

C4(k) 

CgCk) 


.  r)*l  -  k  (1 .  1  -  gi.  , 

C2(k)-|yjs, 

“  s)^+(2k+fl?  r)^]  +  a  ^1  -  I  tan*^  )}  ' 

-  C.Ck)  -  — ^ 

2VS  5  2iBi 

.  jL^f *  in  [(14^  s)2+(2k+fl?  rf]  -  rf  1  -  §  tan*^  , 

SlKtV.  V  1-HfPs  0 

-  tk(lik^'f3)r  +  (2k^.l).], 

.  ^  (kr-.). 


and  finally 


h- ’’A-  °A ■ 


Xg  ■  CgQj^  + 

*3  ■  -  ®6«2  -  V  “* 

X,.2(2k*.l)/^). 


I 
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By  using  these  texns  one  ney  nov  write 


a.  Equation  (3*^)  m 


a  -I- 1^-  1^  -  ^  H. 


(E-l) 


b.  Equations (3*  38)  m 


-  -a  +  IJ  +  |c^  -  Ip  Cg 
and 

c.  the  cooftiinatlon  of  equations  (3*37)  snd.  (3*?0)  as 

*1.^  ■  “>i\  *  Va  i  (1  *  I)  ^  V3  r  *  I)  -  “ia®!. 

*  V5 1  *  1)"^  -  V6  IT  *  i)'^  ■  s  ^ 


Thus, 


,  -X^+X,  Xj 

■  r  it(t/2)  *  i+(s/4)  ■  “  *1  * 


Now,  solving  equation  (E«l)  for  (l-»-2y2)  produces 

k♦(c^,)H 

^*5-— E5—  • 


(E-2) 


(E-3) 


The  ccaiblnatlon  of  (E>2)  and  (E-3)  In  order  to  ellatnate  (l-fZ/2)  yields 
j  ^  [-X  jjk+XjCk-a)  +  HXg  +  OKj^H]  -  aCXj^k+Xg)  -  0.  (E-4) 
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The  solution  of  the  quadratic  equation  Is 


i  1  /r  ’“U 
tr-aj^L - 5~^ 


+  OK, 


/r-kXi.+(k-a)X,  tT  i«ra 

iVM4 - ^  *  -r  >• 

J  (E-5) 


When  (^/U^)  >  0  (X|^  >  O),  the  positive  square  root  must  be  chosen 
because  othexvlse  (^/U^)  approaches  a  finite  limit  as  6  •>  o''^;  this 
Is  of  course  Impossible.  For  (^/U^)  <  0  (Xj^  <  O),  the  negative  square 
root  Is  chosen. 
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APPENDIX  F 
Computations 


The  computations  for  the  present  work  were  accomplished  at  the 
Stanford  University  Con5>utation  Center  on  a  Burroughs  220  Electronic 
Data  Processing  System.  The  installation  at  Stemford  consists  of  a 
Burroughs  Algebraic  Compiler,  the  220  digital  computer  with  10,000 
words  of  core  storage,  and  peripheral  equipment.  The  compiler  accepts 
symbolic  programs  written  in  BALGOL,  an  algorithmic  language  based 
on  ALGOL  [23]  and  produces  machine-language  programs  for  the  computer 
[24]. 

The  programs  listed  in  this  appendix  are  written  in  the  Stanford 
version  of  BALGOL.  The  BALGOL  language  statements  are  reasonably 
self-evident.  In  the  programs  shown,  use  has  been  made  of  the  conqputer's 
ability  to  perform  repetitive  and  complex  arithmetic  calculations  with 
great  speed.  Note  that  each  FOR  statement  indicates  a  group  of 
calcvilatlons  (delineated  by  the  BEGIN  -  END  pair)  which  are  to  be 
repetitively  done.  These  statements  sure  sometimes  nested  within  each 
other.  For  example,  if  a  FOB  group  is  to  be  repeated  10  times  and 
it  contains  another  FOR  group  which  is  to  be  repeated  10  times,  100 
repetitions  will  occ\ir  within  the  second  or  nested  group.  The 
PROCEDURE  SIMPSON  l()  is  a  closed,  independent  routine  which  performs 
numerical  integration  according  to  Simpson's  Rule;  the  procedure  may 
be  called  and  used  at  any  time  in  a  program  in  which  the  procedure 
is  defined. 

When  the  results  are  truncated  by  the  computer  during  calculations 
and  print-out,  no  rounding  occvirs.  Thus,  truncation  causes  an  un¬ 
certainty  of  one  unit  in  the  last  figure  of  all  the  tabulated  data. 
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KpmsDTX.  a 

trcfcvltatlng  Flov  about  a  Slender  Wedge 


Here  the  results  of  Section  3>1  are  used  tw  calculate  the  flov 
parameters  for  a  supercavltating  flov  about  a  slender  vedge  in  a  trans 
verse  gravity  field.  The  correspondlz^  flov  about  a  supercavltating, 
flat- plate  hydrofoil  has  been  studied  by  Parkin  [4],  vhlle  Acosta  [^1 
considered  the  flov  about  a  slender  vedge  in  a  longitudinal  gravity 
field.  The  piirpose  of  this  vork  Is  to  coaqplete  the  studies  of 
gravity  effects  In  linearised  flov.  As  In  Parkin's  vork,  this  treat¬ 
ment  of  the  problem  la  not  capable  of  describing  cavity  flows  vlth 
large,  bouyant  effects.  The  theory  Is  expected  to  be  valid  vhen  the 
effects  of  gravity  are  of  first  order  smallness  consistent  vlth  the 
linearisation  approximations.  The  notation  used  here  Is  consistent 
vlth  that  \ised  previously. 

The  base  flov  Is  an  Irrotatlonal,  Invlseld,  and  Incoavresslble 
\mlform  flov  extending  to  Infinity.  The  upstream  velocity  far  from 
the  vedge  Is  U^,  and  the  origin  of  the  coordinates  (x,y)  Is  taken  at 
the  nose  of  the  unit  length  vedge  vhlch  Is  aligned  syaBmtrlcally  In 
the  flov.  The  reference  elevation  at  infinity  Is,  then,  sero.  The 
flov  Is  sketched  In  Figure  31.  The  acceleration  g  due  to  gravity 
Is  directed  dovnvard  In  the  minus  y-dlrectlon,  perpendicular  to  the 
freestream  velocity  and  vedge  path. 

In  this  flov,  Bemoialll's  equation  is  given  by 

“  Pq  P®^o  “  ^*c  *  *  ^c' 

vlth  the  subscripts  referring  to  Infinity,  the  vedge  surface,  and  the 
cavity  surface  respectively.  Since  the  condition  of  constant  jnressxire 
in  the  cavity  must  be  satisfied  in  steady  flov,  the  non-dimensional 
cavitation  nuidber  may  be  defined  as 

P  “P 
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Another  characteristic  non-dimensional  parameter  peculiar  to  gravity 
flows  Is  the  Froude  Number  F  which  is  given  by 

when  based  on  the  \uilt  length  of  the  wedge. 

La  accordance  with  the  assumptions  of  the  linearized  theory 
developed  In  Section  2,  the  flow  velocities  near  the  body  are  assumed 
to  be  represented  by  the  perturbation  coi^ponents  u  and  v  such  that 


(U  +  u^v'J 


at  any  point  In  the  flow.  Fran  equation  (G-l),  one  has,  to  the  first 
order,  on  the  cavity 


u 

T 


(G-3) 


La  accordance  with  the  basic  assimptlon  of  slenderness  In  the  linearised 
theory,  one  can  argue  that  the  variations  In  the  gravity  term  In 
equation  (G-3)  are  small  over  most  of  the  cavity  (see  Section  2.1). 

Thvis,  the  term  gy  /U  may  be  replaced  by  an  average  term  +  g.  The 
boundary  condition  (G-3)  becomes 


and 


u  o  g 
U  “  2  ■  U  ' 

CO  *0 


(G-3«) 


on  the  cavity  surfaces.  On  the  wedge,  the  boundary  condition  Is  given 
by 


▼(x/Xq) 


“  "  u^+uU^y.)  • 


(0-4) 


If  one  lets  ^1  ~  j  -  a  result  from  the  zero  gravity  case  - 

and  e:]qpands  equation  (G-4)  In  terms  of  U  ,  he  finds,  to  the  first 
order,  that 


and 


(G"Ha) 


V  ■>  OW  ,  y  >  0, 
c 

V  «  -o<U  ,  y  <  0 

C 

on  the  ved^. 

From  the  results  of  Section  2.2,  It  foUovs  that  the  boundary 
conditions  (G-3a)  and  (G-Ua),  together  vlth  three  assusgptlons  -  cavity 
closure,  the  vanishing  of  the  perturbation  velocities  at  Infinity, 
and  smooth  separation  at  the  trailing  edges  of  the  wedge  -  are  sufficient 
to  determine  a  solution  to  a  boundary  value  problem  for  the  complex 
pertiirbatlon  velocity  v.  Recall  that  v  was  defined  In  Section  2.2  as 


w  =  u  -  Iv, 


vlth  V  analytic  outside  the  silt  x-axls  of  the  physical  z-plane. 
The  conqplete  boundary  value  problem  Is  as  follows: 

To  find  v(k.),  analytic  off  the  silt,  such  that 


D.  E. 


B.  C. 


a) 

b) 


c) 

d) 


Real  (v) 
Real  (w) 
Im  (w)  » 
La  (v)  ■ 


dzdz 

V  - 

■ 

1  <  X  <  i,  y  » 

U  a 

=  ”2”  +  8> 

1  <  X  <  i,  y  «  0’ 

-oW^,  05*5^#  y  * 

OU^,  0  <  3t  <  y  -  0’ 


e)  the  cavity  closes,  l.e.,  the  net  source  strength  Is 
zero  on  the  silt. 


f)  w(*)  ^0  as  z  -*-m. 

g)  there  are  no  singularities  at  the  trailing  edges  of 
the  wedge. 
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As  Isefore,  z  b  x  +  ly,  and  the  boundary  conditions  are  applied  on  the 
slit  x-axls  in  the  cosqplex  z»plane. 

If  g  B  -e,  the  above  boundary  value  problem  for  a  "gravity  flow" 
Is  precisely  the  sane  as  that  given  In  Table  1  for  a  uniform  shear 
flow  peist  a  wedge.  Again,  one  sees  the  similarity,  at  least  In  the 
linearized  case,  of  rotational  and  gravity  flows.  Because  of  this 
similarity,  the  results  of  Section  3*1  ney  'be  used  directly  after 
appropriate  changes  of  notation. 

One  finds  first  that  the  gravity  field  has  no  effect  on  the 
cavitation  number  -  cavity  length  equation  or  the  cavity  eurea.  Thus, 

(=-5) 

and  ' 

A^Ba(i3/2.i).  (G-6) 

Second,  on  the  cavity  surfaces  the  horizontal  conponent  of  the 
velocity  Is 


U_(l*f)7i 

on  the  upper  and  lower  surfaces  respectively.  Thus,  one  has  for  the 
cavity  shape 


on  the  \Q)per  surface  and 


on  the  lower  surface.  IVom  equations  (3*13)  of  Section  following 
substitution  of  -g  for  e,  one  finds  that 


I 
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a. 


1K3U 


tan 


t  >  1,  (0-7) 


on  the  upper  surface. 


h. 


|t| 


koSJ 


tan 


oWi?^ 


t  <  -1,  (G-8) 


on  the  lover  surface. 


c. 


1 


L  (t®+T)(t^+R)J 


AlsOj  from  Section  3<^>  the  constants  A  and  D  are 

OSJ 

^  'irrfciT 

and 

‘rS  a.  (vr*  vr). 


(0-9) 


In  the  above, 
and 


T  -  2i  -  1  + 

R  »  2i  -  1  -  211(1-1). 


The  pressure  force  coefficients  are  calculated  next.  The  major 
difference  between  the  rotational  and  gravity  cases  occurs  in  these 
calculations.  !nil8  difference  Is  due  to  the  fact  that  idiUe  tlM  ro¬ 
tational  effect  enters  the  pressxire  coefficient  only  throuc^  the 
perturbation  velocity  terms,  the  gravity  effect  enters  through  both 
the  velocities  and  the  term  pgy  in  the  Bernoulli  equation  (O-l). 

The  pressure  coefficient  C  Is  defined  to  be 

P 
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(G-10) 


C 


P 


Thus, 


and,  from  equation  (G-l),  one  obtains 


C  =  CT  +  1 
P 


lo  20^0 

T'T- 


In  accordance  with  Parkin's  dlacusslon[^]  regarding  the  linear  contrlbu- 

2 

tlon  from  the  second  order  term  (u/U  )  ,  one  obtains  after  linearization 

00 

the  relation 


Where  +  replaces  sy^/^  cn  the  upper  and  lover  wedge  surfaces 

respectively.  One  should  refer  to  Section  2.3  for  a  detailed  discussion 
of  this  substitution  and  should  conq;>a3re  equation  (G-U)  with  Its  rota¬ 
tional  equivalent  equation  (3.IU).  By  using  the  results  of  that  section, 
one  has 

-  i-  .  1  +  in  (yr+  fPDtaaeJ^  (O-Ua) 

Q  p 

With  X  s  i  cos  d/(|-sln  d).  The  minus  sign  and  0  <  6  <  TLf2  apply  to 
the  upper  wedge  surface;  the  plus  sign  and  II/2  <  d  <  ^  apply  to  the 
lower  surface. 

The  remainder  of  the  force  coefficients  are  determined  from 
equation  (O-U).  First,  the  drag  coefficient  is  given  in  Table  4 
as 


D 

ipuf(2a) 


(G-12) 


BODY 
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Since  U  dy  ■  vdx  on  the  wedge  and  2uv  »  -  Ik  w^(z),  one  may  write 
0 


-  dz  +  2^+  xdy  +  (2+a)|^dy  . 

JBODY  »  BODY  '  *  BODY 


From  Section  3*1 *2^ 


dz  «  -  Im  J, 


BODY  « 


,j,  +  (2+a)^+  dy  +  (2+(r)  f^dy- 


The  closure  condition  for  this  cavity  model  requires  that 


By  using  these  res^llts,  one  obtains 


b  ^dy  e  0. 
BODY  +  CAV 


since  X  B  y/a  on  the  wedge.  Finally,  it  is  easily  shown  that  the 
last  two  contour  Integrals  in  the  above  are  zero;  hence,  with  the 
value  of  given  in  Appendix  B, 


(2+cr)^a  i 

. •  n 


(0-13) 


Second,  the  lift  coefficient  is  defined  in  Table  U  as 
Cj  =  — -  B®C  dx. 

iP«f(CS0B«B0OT  ^ 

From  Equation  (G-U), 


'1.-2 


After  one  performs  the  indicated  integrations  and  notes  that 
Re(w/U^)dz  »  (u/U^)dx  on  the  slit,  the  above  equation  becomes 


Cj^  ■  2^^^  -  (2+0)  Re^^  dz. 
'  *'  BODY* 
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The  results  of  Section  3*1 *2  show  that 


Re^3^  dz  »  ^  in  (VT  +  -2(i-l)l  . 

body  *  "  ^ 


Thus, 


^  ^i-(2+(T)^Vi(i-iy  ji  +  in(i(r+  ifliiyj  -  . 


(0-14) 

Third,  the  moment  coefficient  about  the  nose  of  the  ved^  Is 
defined  In  Table  4  as 


.  L-(dl.t.  to  L) 

•*“  iaficmsD)^ 


To  the  first  order. 


=  2  J'  x^dx  -  J'  x^dx  I  -  (2+ct)  ^  u  xdx. 

“\o  1  /  BODY 


Following  simplification  and  introduction  of  complex  notation,  this 
equation  becomes 


°i«>  ■ 

'  *'  BODY 


In  Section  3*1«2,  It  Is  shown  that 


Re^wzdz  -  ^  j\8{i-lj(2i+l)  ^ 

.2(i^-l)^  . 


Thus, 


4Tin(yiiyC?)/  2fi(i-iy(l4Tl 


(!r-l)(l-T) 

+  (1-T)  - 


-  2(i^-l)  J  . 


(0-15) 
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Filially,  It  remains  to  define  the  gravitation  parameter  g/U  .  On 

00 

the  hasis  of  the  analogy  between  the  linearized  rotational  and  gravity 
problems,  this  definition  is  taken  directly  from  equation  (3.23)  of 
Section  3«l«2j  hence, 

t  '  (§)  51^  • 

By  introducing  equation  (G-6)  and  the  Froude  Nusiber  F  into  this  result, 
one  has 

00 

The  complete  solution  is  suneiarlzed  in  Table  9* 

Certain  results  of  the  gravity  flow  analysis  are  the  same  as 
those  previously  obtained  in  the  study  of  uniform  shear  flow  past  a 
wedge.  These  duplicated  results  incl\]de  the  cavity  length- cavitation 
n\miber  relationship,  the  cavity  area,  and  the  drag  coefficient  -  ell 
of  which  are  Independent  of  gravity  effects.  Figures  7  and  6  show 
the  drag  coefficient  and  cavity  area  as  functions  of  cr.  The 
i-o  relationship  is  plotted  in  Figures  I6  and  17;  the  curve  for  gravity 
flow  is  the  same  as  that  for  Irrotatlonal  gravity -free  flow,  i.e.,  the 
curve  labelled  e/u^  s  0.  As  in  the  gravltywfree  case,  the  cavity  length 
is  limited  in  this  linearized  theory  so 


-  1) 

2^(i-l) 


(G-16) 


It  is  possible,  however,  that  as  i  •  (a  O),  the  effects  of  gravity 
may  be  exaggerated.  Finally,  since  i/<  is  eq:ual  in  magnitude  to 
Figure  9  gives  a  plot  of  both  parameters. 

The  remainder  of  the  numerical  results  are  listed  in  Appendix  F 
and  Illustrated  in  Figures  32  through  37  •  Figures  32  through  3? 
picture  the  pressure  force  coefficients.  The  first  two  figures  show 


the  effect  of  cavitation  nuaiber  on  and  in  a  gravity 

Note  that  both  coefficients  are  linear  functions  of  l/f^.  Hence, 


independent  of  a,  both  coefficients  Kre  inverse  funetlcms  of  the 
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square  of  the  Frouda  nusiber.  As  a  approaches  zero,  the  lift  and 
moment  become  negatively  Infinite  for  all  finite  Froude  numbers.  The 
limitations  of  the  theory  are  seen  clearly  In  Figures  3^  uid  35^  which 
show  the  pressure  coefficient.  When  the  Froude  nuniber  Is  small,  the 
gravity  effects  are  large  and  the  pressure  coefficient  on  the  lower 
side  of  the  wedge  Is  negative  over  the  whole  body  surface.  Such  a 
condition  (one  In  which  the  pressure  on  the  body  Is  always  less  than 
the  pressure  In  the  cavity)  Is  contradictory.  However,  as  In  the  case 
of  uniform  shear  flow  past  a  wedge.  It  seems  permissible  to  allow 
a  negative  pressiire  In  the  Immediate  vicinity  of  the  nose  of  the  wedge 
for  the  reasons  presented  In  Section  3*1*3«  Thus,  In  Figure  35 
for  excuDcle,  the  curve  of  C  vs  x  for  F^  s  i6  Is  a  reasonable 
approximation,  \dille  when  F^  =  4  It  Is  seen  that  the  lladts  of  the 
theory  have  been  exceeded.  This  behavior  corresponds  to  that  found  by 
Parkin  [4]  In  his  study  of  gravity  effects  on  hydrofoils.  Li  the 
present  case,  the  behavior  Is  due  to  the  Increased  size  of  the  gravity 
effects,  represented  by  g/U^,  for  small  Froude  nusdbers.  Since  this 
theory  was  expected  to  hold  only  for  flows  with  small  gravity  effects. 
Its  use  must  be  restricted  accordingly.  The  behavior  of  acts, 
then,  as  a  guide  to  the  limits  of  the  theory. 

The  final  figures.  Figures  36  and  37f  show  the  effect  of  Froude 
number  on  cavity  shape  and  the  effect  of  cavitation  nuidber  on  the 
location  of  the  center  of  lift  at  an  arbitrary  FToude  nuid^er  (since 
Cj^/Cl  Is  not  a  function  of  F).  The  cavity  shape  Is  seen  to  be 
distorted  downward  In  the  middle  and  upward  at  the  end.  The  effect 
of  the  transverse  gravity  field  Is  exactly  opposite  to  the  effect  of 
a  uniform  shear  flow  with  positive  vortlclty  (see  page  33  and  Figure 
10).  The  fact  that  the  cavity  Is  not  inclined  upward  by  bouysney  Is 
no  longer  surprising  since  the  same  result  was  predicted  by  Parkin'  s 
analysis  and  has  been  confirmed  by  e^^rlments  on  cavltatlmi  behind 
tWD-dlmensloiMd  bluff  bodies  at  the  California  Institute  of  Technology 
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d.  U(«)  -  ortJ^,.  0  <  K  <  1,  d.  !•(•)  -  0  <  *  <  I, 

y  -•  O’  y  -•  0* 

«.  Th«  cavity  it  elottd,  i.«. ,  the  net  tenrce  strength  on  the  elit  is  tero. 

f.  •(«)  0  at  t  i.e.,  (u,v)  ■  0  at  infinity. 

g.  «(i)  nuat  not  contain  nonintegreble  aingnlaritiea  on  the  alit  or  have 
■altiple  values  off  the  alit. 

h.  The  flow  ia  characterised  by  a  anooth  aaparation  fron  tha  rear  of  tha  body, 

i.  e. ,  w  <  •  at  X  ■  1,  y  -•  0. 
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!■(*)  .  ceastaat  e.  foil 


TABLE  4 


PRESSURE  FORCE  COEFFICIENTS 


Coefficient 


P  •  P 

1.  C.  -  - ^ 


Wedge 


(2+2)  JL  *  +  5  .  JL 


V  U„  2 


U„/ 


Hydrofoil 

.(2+2)f  X  X  +  -S-  . 

V  U„  U«  2  / 
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Besed  on  firet-order  emellnees  of  engle  of  eCteck  ou 
**Teke*  et  lending  edge  of  body;  poaitire  in  tke  connter-clocknise  direction. 
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TABLE 


The  •■bscripts  U  ud  L  r«f«r  r«*p«ctiT«ly  to  tke  oppor  ood  looor  aorfaccs  of  tke  body  or  cavity. 


TABLE  6.  SUMMABY  OP  RESULTS  FOR  ASYMMETRIC  WEDGE  FLOW. 


OuMtity 

EqvatioB 

No. 

Cavity 

Laiifth 

2:  .  a  /  ^  v5_M  +  2^3  \ 

A  -  1  ^  J 

(3.8) 

Cavity 

Araa 

A^  -  -  1) 

(3.11b) 

Cavity 

Shape 

(See  equatioBB  3.13  a,  b,  aad  c) 

(3.13) 

Preaaure 

CeafficieBt 
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L’^  \  l-l  1  ■  aie  e ) 

Upper 

Serface 

+  -  1  +iLkII  tee  dAm  (/? +/frr)  +  x^l 

U. \w  7r(l-»T)  /J 

(3.16) 

0  <6  <n/2 

Lover 

Serfaee 

+  ^  •  1  +  tae  6  •  'tm  (A+A^)  •  x^l 

U,^w  ffd+T)  JJ 

(3.16) 

7r/2  <6<v 

Drag 

CoefficioBt 

Lift 

CoefficioBt 

,  _  -t  COB^  6 

't-  aiB*  9 

Cjj  ■  (2  +2)2  orC/wl-C  -  1 ) 

(3.15) 

(3.18) 

Cl  J(2+2)^/l(^i)  1^1+^  4  {V^+v^)J  +i 

CjK,  -(xV2+Z)proT(li^) 

(3.20) 

Moaeet 

CoefficioBt 

.  T  xa/tf-t-i)  (T+1)1\ 

'  (t2.1)(1.T)  L  J/  3 

--^2^  (3.22) 

Vorticity 

Paraaieter 

X  -  /  € 

Uo)  Vu./ 

(3.23) 

NOTE:  T  - 

1  +2  /ti-t-i). 
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TABLE  7. 


SUMMARY  OF  RESULTS  FOR  SYMMETRIC  WEDGE  FLOW 


Quantity 

Equation 

No. 

Cavity 

Length 

1  fa  +f  M  .  a  All  +  lA  \ 

2+crL  \U,  /  l-l  1  v\  1 

(3.26) 

Cavity 

Area 

A^  =  -  1) 

(3.28b) 

Cavi ty 

Shape 

..aJfe/fc.i.  4(^-1)  •] -JlillULlL,  4>l 

W  J  L  c  CJ  (^2+x)2({2+n)2 

(3.29a) 

X  -  Hi  .4('t-l)t2/  (t^+TXt^+R)] 

(3.13c) 

Preaaure 

Coefficient 

C,  .  ,2  4.((X)(..,  ||*-«|)] 

0  <  0  <  w/2 

(3.32) 

X  ■  't  coa^  6  /  (“1  •  ain^  6) 

(3.31) 

Drag 

Coefficient 

Cjj  -  (2  +  )*  at/7r('t-l) 

(3.33) 

Vorticity 

X  .(  € 

U,  VU,/  2('t-l) 

(3.  34) 

NOTE:  T  - 

2t-l  +  2/f(^l):  R  “  a^-l-2 

I 
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TABLE  9. 

Quantity 

Cavity 

Length 

Cavity 

Area 

Cavity 

Shape 

Presaure 

Coefficient 

Upper 

Surface 


Lower 

Surface 


Drag 

Coefficient 

Lift 

Coefficient 


Moaent 

Coefficient 


Gravi ty 

Paraneter 


NOTE:  T  - 


SUMMARY  OF  RESULTS  FOR  lEDGE  FLOW  IN  TRANSVERSE  GRAVITY  FIELD 


Equation 

No. 

• 

_OL  =  a  vi+1  +  2J^  ) 

2+cr  7r\  /J.l  't-l/ 

(G-S) 

- 

A^  ■  -  1 ) 

{G.6) 

(See  equationa  G>7,  8,  and  9) 

(G.7,8,9) 

■  -'(it)  “  *  'r:i: «') 

.  .  1  /Jn  (  /r  +  /t-l )  tan  ^  \ 

U^Itt  yrd-PT) 

(6-lla) 

0  <  ^  <  7T/2 

Cp,  .  «( a) .  *  ,JW{2  +4  |»  •-  e|) 

A 

.  Jl\M  .  1  +  Midi  ‘Ik  (v^+/CT)  tan 

U.Lw  Trim) 

* 

Tr/2  <  ^  5  w 

• 

_  -fc  coa‘  6 

X-ain^  9 

Cq  •  (2-K^)^  artAr(^-l) 

(G-13) 

«  2 -JL  -(2-fO')  +  — --'Cn  +/t- 1  )j  -'t+1^ 

(G.U) 

c^.jil i.iiisifmT)  (»1«)  fi  ♦ 

^  0«(_3  2  ^  L  (T*-1)(1.T) 

(G-IS) 

.  ^  (1.4T)  +  i.xjj  -2(^*.l))j 

-• 

i  -  o(<t3/2.l)/2F2('t.l) 

(G-16) 

ai-i  +  2  VlfJLx) 


I 
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UNIFORM  FLOW 

U„ 


»  FERTORBED  FLOW 


a-  Flova 


I  WEBOES 


2.  HYDROFOIL 

b.  Probleaa 


FIGURE  1.  FLOW  PATTERNS  AND  PROBLEMS. 


b.  Hy4rofQil 


FIGURE  2.  FULLY  CAVITATED  FLOWS. 


b.  Hydrofoil  flow 

FIGURE  3.  MAPPING  OF  t-PLANE  ONTO  UNIT  CIRCLE. 


I 
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FIGURE  4.  SYMMETRIC  SHEAR  FLOW  PAST  A  WEDGE. 


C-PUNE 

SyMctric  wedge  flow 


FIGURE  S.  MAPPING  OF  s-PLANE  ONTO  UNIT  CIRCLE. 


FIGURE  1.  CAVITY  AREA  VS  CAVITATION  NUMBER  FOR  FLOW  PAST  A 
UNIT  WEDGE. 


-  112  - 


CA,¥lTATION  NUMBER  >  £  (OR  a) 


FIGURE  9.  VORTICITY  PARAMETER  OR  GRAVITY  PARAMETER  VS 
CAVITATION  NUMBER  FOR  FLOW  PAST  A  UNIT  WEDGE. 
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riOURC  10.  CAVITY  SHAPES  AT  CONSTANT  CAVITATION  NUMBER  IN  UNIFORM  SHEAR  FLOW  PAST  A  UNIT  WEDGE 


100  0 


(•n/3)/''3  -  1II3IDIJJ30D  an 


o 

h. 


H 


i 

i 


-  115  - 


FIGURE  11.  LIFT  COEFFICIOIT  FOR  SHEAR  FLOW  FIGURE  12.  HOMENT  COEFFICIENT  FOR  SHEAR 

PAST  A  WEDGE.  PAST  A  WEDGE. 


raESSVRE  COEFFICIENT  -  C.  PRESSURE  COEFFICIENT 


0  0.2  C.X  0.6  0.8  1.0 


DISTANCE  FROM  NOSE  -  x 

FIGURE  13.  1  ”  0.086S,  ^  =  30,  a  -  10®  PRESSURE  COEFFICIENT 

FOR  SHEAR  FLOW  PAST  A  WEDGE. 


FIGURE  14.  I  ^  0.0483,  'C  90.  a  -  10®  PRESSURE  COEFFICIENT 
FOR  SHEAR  FLOW  PAST  A  WEDGE. 


I 
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Noimuvo 


CAVITATION  NUMEI  -  2 

FIGURE  15  CENTER  OF  LIFT  FOR  SHEAR  FLOW  PAST  A  VEOOE. 


I  2  AOIIO  20  NOIOOOlOO  200  NOOOOO 


CAVITY  LENOTN  -  'C 

FIGURE  16.  CAVITATION  NUMBER  VS  CAVITY  LENGTH  FOR  SYMMETRIC  SHEAR 
FLOW  PAST  A  WEDGE  o  •=  S“. 
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’  0-2 
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2  0.08 

t  0  06 


FIGURE  17. 
FLOW  P 


l/U, -0.08^  O.M^|e.(^ 

6  8  lO  20  80  80  80  100  200  800  600  iOOO 

CAVITY  LENOTN  - 't  •O® 

CAVITATION  NUMBER  VS  CAVITY  LENGTH  FOR  SYMMETRIC  SHEAR 
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i  “ ' 


ui  0.8 

8 


rii 


0  O.i  0-2  0.3  0.8  0.5  0,8  0.7  0.8  0.8  1.0 

CAVITATION  NUMBER  -<r 

FIGURE  18  DRAG  COEFFICIENT  VS  CAVITATION  NUMBER  FOR  SYMMETRIC 
SHEAR  FLOW  PAST  A  WEDGE. 
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I 


O.OM  0.006  0.008  0.010  0.02 

CAViTATlOi  NUM8ER  -  1 

FIGURE  20.  CAVITY  LENGTH  VS  CAVITATION  NUMBER  FOR 
SHEAR  FLOW  PAST  A  HYDROFOIL:  a  >  1°. 
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CAVITY  LQI6TH 


0.01  0,02  0,00  0,00  O.W  0,1  0.2  O.S 

CAVITATION  NiIMKR  -  £ 


FIGURE  21.  CAVITY  LENGTH  VS  CAVITATION  NUMBER  FOR  SHEAR 
FLOW  PAST  A  HYDROFOIL:  a  -  4**. 
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FIGURE  22.  EFFECT  OF  VORTICITY  AND  CAVITATION  NUMBER  ON 
VORTICITY  PARAMETER  FOR  SHEAR  FLOW  PAST  A  HYDROFOIL. 


FIGURE^  23  CAVITY  SHAPES  AT  CONSTANT  LENGTH  IN  UNIFORM  SHEAR  FLOW. 


NORMAL  FORCE  COEFFICIENT 


0  0.1  0.2  0.8  0.0  0.5 

CAVITATION  NUMOER  -  Z 


FIGURE  24.  NORMAL  FORCE  COEFFICIENT  VS  CAVITATION  NUMBER 
FOR  SHEAR  FLOW  PAST  A  HYDROFOIL. 
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DISTANCE  TO  CENTER  OF  LIFT  -  x  NONENT  COEFFICIENT 
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0.1 
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-0.02  0 

0.02 
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0-08 

RELATIVE  VORTICITY  - 

f/0« 

FIGURE  2' 
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FIGURE  28  EFFECT  OF  VORTICITY  ON  *  AT 
VARYING  I  FOR  FLOW  PAST  A  HYDROFOIL 
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FI^RE  32  EFFECT  OF  CAVITATION  NUMBER  ON  LIFT  FIGURE  33.  EFFECT  OF  CAVITATION  NUMBER  ON  MOMENT 

COEFFICIENT  IN  GRAVITY  FIELD.  COEFFICIENT  IN  GRAVITY  FIELD 


ntESSUIE  COEFFICIEMT  -  C.  PRESSURE  COEFFICIENT 


FIGURE  35.  CT  ^  0  0483,  I  -  90,  o  =■  10°  EFFECT  OF  FROUDE 
NUMBER  ON  FOR  FLOW  PAST  A  WEDGE. 
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iiVNIOM  AilAVa 


FIGURE  37.  EFFECT  OF  CAVITATION  NUMBER  ON  ¥  FOR  ARBITRARY  FROUDE 
NUMBER  FOR  FLOW  PAST  A  WEDGE. 
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